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Pre face 


The  study  of  aircraft  flight  control  is  a  very  broad  field.  In 
selecting  a  thesis  topic,  I  was  interested  in  choosing  a  project  that 
would  permit  me  to  experience  as  many  aspects  of  this  field  as  possible 
This  necessitated  that  the  project  be  small  enough  in  scope  that  it 
afforded  me  the  opportunity  to  engage  in  analysis  as  well  as  design. 

The  NEXTRPV  program  was  well  suited  to  meet  this  goal.  Whereas  this 
project  is  not  necessarily  glamorous,  it  gave  me  the  chance  to  exper¬ 
ience  "basic"  engineering  which  I  believe  is  on  important  part  of  any 
education. 

This  thesis  originally  began  as  strictly  a  flight  control  system 
analysis  and  design  for  an  existing  remotely  piloted  vehicle  (RPV).  As 
the  project  progressed,  it  became  evident  that  the  RPV  was  fairly  stabl 
and  would  not  need  as  much  flight  control  augmentation  as  initially 
anticipated.  Therefore,  this  thesis  was  extended  to  include  an  investi 
gation  of  a  practical  application  of  this  RPV. 

Many  people  contributed  to  this  thesis.  To  each  of  them,  I  say 
thank  you.  Special  thanks  goes  to  my  thesis  advisor,  Captain  James  T. 
Silverthorn,  for  his  numerous  suggestions  and  neverending  attention. 

My  gratitude  is  extended  to  Professors  Robert  A.  Calico  and  John  J. 
D'Azzo  whose  assistance  was  invaluable.  My  appreciation  is  also  ex¬ 
pressed  to  Mrs  Anna  Lloyd  for  preparing  this  manuscript.  Lastly,  I 
would  like  to  thank  my  wife,  Susan,  whose  constant  support,  encourage¬ 
ment,  and  understanding  made  the  completion  of  this  project  possible. 

Brian  L.  Jones 
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Abstract 


~  “This  study  investigated  the  stability  and  control  of  a  remotely 
piloted  vehicle  (RPV)  with  a  thrust  vectoring  unit  attached.  All 
geometric  and  aerodynamic  data  was  generated  and  used  to  analyze 
the  RPV.  Specific  handling  qualities  were  developed  and  compared 
with  the  RPV  characteristics.  This  comparison  indicated  that  the  RPV 
was  too  oscillatory  in  both  the  phugoid  and  dutch  roll  modes.  Also, 
the  RPV  displayed  a  dominant  spiral  mode.  A  flight  control  system 
was  synthesized  to  eliminate  these  traits.  Evaluation  of  this  flight 
control  system  was  conducted  through  the  use  of  three  different 
sensitivity  studies  and  a  nonlinear  simulation.  In  addition,  a  model 
matching  application  was  examined  for  this  RPV.  Model  matching  entails 
using  design  procedures  to  synthesize  an  expanded  flight  control  system 

so  that  the  RPV  has  dynamic  characteristics  similar  to  the  F-15. 

K 
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FLIGHT  CONTROL  SYSTEM  ANALYSIS  AND  DESIGN  FOR 
A  REMOTELY  PILOTED  VEHICLE  WITH  THRUST  VECTORING  UNIT 

I.  Inf.  reduction 


Background 

In  aircraft  design,  the  aeronautical  engineer  is  faced  with  many 
design  objectives.  For  military  aircraft,  these  objectives  would  in¬ 
clude  such  traits  as:  low  observability,  high  lethality  and  survivabil¬ 
ity,  STOL  capability,  and  lastly,  air-to-air  and  air-to-ground  capabil¬ 
ity.  Several  schools  of  thought  exist  on  how  to  satisfy  these  objectives. 
One  alternative  is  the  use  of  a  thrust  vectoring  unit  (TVU).  With  such 
a  unit  employed,  the  tail  section  could  be  removed  without  any  reduction 
in  aircraft  control .  This  would  reduce  the  infrared  signature  and  hence 
lower  an  aircraft's  observability.  Also,  thrust  vectoring  increases  an 
aircraft's  maneuverability,  which,  in  turn,  increases  its  lethality  and 
survivability  while  enhancing  its  air-to-air  and  air-to-ground  caoabil- 
ity.  In  addition,  thrust  vectoring  has  been  shown  to  improve  an  air¬ 
craft's  STOL  capability. 

With  a  continual  reduction  in  funding,  the  Air  Force  needs  to  find 
inexpensive  ways  in  which  to  develop  and  test  such  theories  as  thrust 
vectoring.  One  answer  to  tin's  problem  is  using  remotely  piloted  ve¬ 
hicles  (RPVs)  as  flight  test  vehicles.  RPVs  can  be  operated  at  a  small 
fraction  of  the  cost  of  a  test  aircraft,  resulting  in  substantial 
savings  in  research  dollars. 

The  two  above  mentioned  concepts  were  combined  to  form  the  NEX1RPV 
program  under  the  direction  of  the  Control  Systems  Development  Branch 


(AFWAL/FIGL ) ,  Flight  Dynamics  Laboratory,  Air  Force  Wright  Aeronautical 
Laboratories.  Their  stated  objective  was  to 

demonstrate  the  use  of  practical  thrust  vectoring  and  digital 
flight  control  systems  in  a  low  cost  flight  research  vehicle 
oriented  toward  a  more  stealthy  design  with  an  expanded  maneuver¬ 
ing  envelope. 

The  approach  the  laboratory  intended  to  pursue  was  to  first  integrate 
the  TVU  into  the  RPV  through  an  analog  control  system.  Flight  control 
authori ty  would  gradually  be  transferred  to  the  TVU  over  a  series  of 
test  flights.  Following  the  successful  completion  of  the  TVU  full 
authority  test  flight,  the  entire  procedure  would  be  repeated  using  a 
digital  control  system.  At  the  conclusion  of  this  phase,  the  tail 
section  would  be  removed  for  a  final  performance  evaluation. 

The  RPV  chosen  for  this  program  suited  it  well  for  several  reasons. 
Since  this  RPV  was  obtained  from  a  previous  program,  it  had  already  been 
flight  proven.  Thus,  all  considerations  other  than  flight  control  system 
(FCS)  design  were  eliminated  from  this  project.  Also,  the  twin  boom/ 
vertical  tail  assembly  required  little  modification  for  the  TVU  addition 
(see  Fig  ").  Lastly,  since  the  RPV  wa s  already  constructed,  the  associ¬ 
ated  development  costs  were  saved. 

Problem  Statement ,_Scope_,_and  Assumptions 

The  purpose  of  this  study  was  to  determine  the  stability  and  con¬ 
trol  characteristics,  develop  specific  handling  qualities,  synthesize  a 
FCS  to  satisfy  these  handling  qualities,  and  investigate  a  practical 
application  of  this  RPV. 

This  study  was  restricted  to  the  TVU  zero  authority  case.  In  other 
words,  the  TVU  was  attached  to  th<^  RPV,  but  was  given  no  control 
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authority.  The  only  exception  appeared  in  the  basic  analysis.  In  this 
phase,  center  of  gravity  location  and  TVU  attachment  were  used  as 
parameters.  X-axis  center  of  gravity  (xcg)  assumed  three  locations: 

10,  25,  and  40  percent  mean  aerodynamic  chord  (mac).  The  TVU  attached 
and  unattached  cases  were  examined  for  each  xcg  location.  The  practical 
application  was  restricted  to  the  longitudinal  case.  As  shown  in  Chapter 
V,  the  lateral -directional  case  would  simply  be  a  repetition  of  the 
same  theory. 

Initially,  several  simplifying  assumptions  were  imposed.  Perturba¬ 
tions  about  straight  and  level  flight  were  assumed.  This  allowed  the 
linearized  aircraft  equations  of  motion  to  be  decoupled  into  longitudinal 
and  lateral-directional  sets.  For  the  FCS  design,  only  pitch,  roll,  and 
yaw  rates  were  acceptable  feedbacks.  It  was  further  assumed  that  a 
processing  unit  would  be  accessible  should  compensation  be  required. 
Lastly,  full  state  feedback  was  assumed  available  for  the  practical 
application. 

Approach  and  Presentation 

This  study  was  separated  into  four  major  phases.  Phase  one  con¬ 
sisted  of  obtaining  all  necessary  geometric  and  aerodynamic  data  for 
this  RPV ,  followed  by  the  generation  of  six  sets  of  stability  derivatives 
(TVU  attached  and  unattached  cases  for  three  xcg  locations).  These  sets 
made  up  the  mathematical  models  used  in  this  study.  Chapter  II  of 
this  report  presents  the  results  of  phase  one. 

The  basic  analysis  of  the  RPV  constituted  phase  two.  This  analysis 
began  with  the  development  of  specific  RPV  handling  qualities.  The 
aircraft  equations  of  motion  were  derived  and  a  modal  analysis  of  the 
six  models  generated  in  phase  one  was  conducted.  Phase  two  concluded 
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with  a  comparison  of  the  modal  analysis  with  the  RPV  handling  qualities. 
The  report  on  phase  two  is  contained  in  Chapter  III. 

Phase  three  involved  the  design  and  evaluation  of  a  FCS  to  satisfy 
all  the  handling  qualities.  Sensitivity  studies  and  a  nonlinear  simu¬ 
lation  were  used  to  evaluate  the  FCS.  Phase  three  is  documented  in 
Chapter  IV. 

Lastly,  phase  four  entailed  an  investigation  into  a  practical  appli 
cation  (model  matching)  for  this  RPV.  An  expanded  FCS  was  synthesized 
to  make  this  RPV  have  dynamic  characteristics  similar  to  an  F-15.  Three 
different  design  procedures  were  examined  and  are  described  in  Chapter  V 

The  final  chapter  of  this  report,  Chapter  VI,  presents  conclusions 
and  recommendations. 
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II.  RPV  Modal  Genera  Mon 

Aircraft  stability  and  control  characteristics  are  usually  examined 
based  upon  a  linearized  set  of  differential  equations  of  motion.  The 
development  of  these  equations  (see  Chapter  III)  contains  numerous  param¬ 
eters  which  fall  into  three  categories:  physical  quantities,  flight  con¬ 
dition  defining  quantities,  and  stability  derivatives.  The  purpose  of 
this  chapter  is  to  generate  these  parameters  for  this  particular  RPV. 

Method 

The  physical  quantities  (mass,  moment  of  inertia,  etc.)  and  flight 
condition  defining  terms  (airspeed,  dynamic  pressure,  etc.)  were  easily 
determined  from  the  size  of  the  RPV  and  its  expected  flight  regime. 

Thus,  the  major  crux  of  the  model  generation  became  calculating  the 
stability  derivatives.  Applicable  publications  were  consulted  in  order 
to  determine  the  most  appropriate  method  of  calculating  these  deriva¬ 
tives.  Most  publications  reviewed  listed  The  USAF  Stability  and  Control 
DATCOM  (Ref  8)  as  a  reference.  These  manuals  contain  empirical  formulas 
and  methods  for  determining  stability  and  control  derivatives,  aircraft 
moments  of  inertia,  etc.  A  computerized  version  of  these  methods  appears 
as  Ill®.  USAF  Stabil  i  ty  and  Control  Digi  tal  Datcom  (Ref  18) .  Due  to  the 
widespread  acceptance  of  Datcom  methods  and  the  versatility  of  a  com¬ 
puter  program,  it  was  decided  that  Digital  Datcom  would  be  the  best 
source  of  stability  derivatives. 

Input  Data  for  Digi tal  Datcom 

Digital  Datcom  requires  three  types  of  input  data:  aircraft 
geometry,  airfoil  section  characteristics,  and  flight  condition. 


Aircraft  geometry  was  obtained  by  directly  measuring  the  RPV's  physical 
dimensions.  A  circular  cross-section  approximation  was  employed  for 
the  fuselage  and  booms  (see  Fig  1).  Table  I  is  a  summary  of  all 


Figure  1.  Circular  Cross-Section  Approximation 
for  the  Fuselage  and  Boom 


measurements. 

Two  major  modifications  were  planned  for  the  RPV  and  were  incor¬ 
porated  into  the  mathematical  model.  The  first  was  a  two  foot  nose 
extension  and  the  second  was  the  addition  of  the  TVU.  Both  modifica¬ 
tions  were  modelled  as  part  of  the  fuselage. 

Next,  airfoil  section  character-; sties  were  required  for  each  air¬ 
foil.  The  wing  was  found  to  have  a  varying  airfoil  section,  a  NACA 
23012  at  the  root  and  a  NACA  441,?  at  the  tip.  The  horizontal  tail  was 
a  NACA  0009  while  the  vertical  tails  were  NACA  0010  airfoil  sections. 
Since  NACA  airfoil  sections  were  Found  (Ref  1),  Digital  Dotcom's  in¬ 
ternal  aerodynamics  package  was  utilized  for  the  generation  of  section 
characteristics.  An  approximation  of  the  wing's  varying  airfoil  section 
was  formulated  by  computing  two  sets  of  stability  derivatives,  one  set 
for  each  of  the  root  and  tip  airfoil  sections,  and  linearly  averaging 
the  results.  This  approximation  should  be  accurate  since  the  differ¬ 
ence  between  the  individual  stability  derivatives  was  small. 
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Table  I.  Physical  Dimensions  of  the  RPV 
a.  Wing,  Horizontal  Tail,  and  Vertical  Tail  Dime 


W. 

H.T. 

Apex,  ft: 

x-axis 

4.25 

11.39 

z-axis 

0.60 

2.10 

Chord,  ft: 

root 

2.52 

1 .29 

tip 

2.00 

1 .29 

Incidence  Angle,  deg: 

0.95 

0 

Span,  ft: 

exposed 

13.00 

4.88 

theoretical 

13.71 

5.00 

Sweptback  Angle,  deg: 
leading  edge: 

inboard  panel 

3.50 

0 

outboard  panel 

3.50 

0 

trailing  edge: 

inboard  panel 

0 

0 

outboard  panel 

-1 .40 

0 

*  See  Ref  18  for  definitions  of  terms. 
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Table  I.  (continued) 


b.  Fuselage  Dimensions 

J 

1  x-station,  ft 

perimeter,  ft 

? 

area,  ft 

radi us , 

ft 

2.00 

2.42 

2.83 

3.25 

3.67 

1  3.83 

4.25 
!  4.75 

5.25 

5.75 

6.25 

6.75 

3.33 

3.42 

3.58 

3.71 

3.79 

3.83 

3.75 

3.75 

3.75 

3.75 

3.75 

3.75 

0.88 

0.93 

1.02 

1  .09 

1.14 

1.17 

1.12 

1.12 

1.12 

1.12 

1.12 

1.12 

0.53 

0.54 

0.57 

0.59 

0.60 

0.61 

0.60 

0.60 

0.60 

0.60 

0.60 

0.60 

c .  Boom 

Dimensions 

|  x-station,  ft 

perimeter,  ft 

area,  ft2 

radius, 

ft 

5.00 

0.67 

0.035 

0.11 

5.42 

1  .58 

0.20 

0.25 

5.83 

1  .83 

0.27 

0.29 

6.33 

1.92 

0.29 

0.30 

6.75 

1.67 

0.22 

0.26 

7.17 

1  .62 

0.21 

0.26 

7.58 

1.58 

0.20 

0.25 

8.00 

1.50 

0.18 

0.24 

8.42 

1 .42 

0.16 

0.22 

8.83 

1 .35 

0.15 

0.22 

9.25 

1.29 

0.13 

0.21 

1  9.67 

1  .21 

0.12 

0.19 

1 

I  10.03 

1.12 

0.10 

0.18 

10.50 

1  .00 

0.080 

0.16 

10.92 

0.92 

0.067 

0.15 

! 

11.25 

0.83 

0.055 

0.13 

8 


Flight  condition  defining  parameters  was  tr n  last  sot  of  input  data 
needed.  After  consideration  of  the  engine  capabilities  and  the  flight 
profiles  expected,  the  airspeed  regime  was  determined  to  range  from 
Mach  0.05  to  Mach  0.20.  Similarly,  the  angle  of  attack  range  was 
formulated  as  -4  deg  to  15  deg.  Sea  level  conditions  were  assumed 
for  all  quantities  requiring  atmospheric  conditions. 

Fonviul  it  ion  of  the_  RPV  Representati  ons 

The  choice  of  Digital  Datcom  specified  the  basic  RPV  representation 
as  a  straight  tapered  wing  and  one  vertical  tail  of  increased  area.  It 
was  decided  to  employ  a  scale  factor  of  1.75  in  increasing  the  vertical 
tail  area.  This  choice  was  somewhat  arbitrary  since  no  documentation 
could  be  found  on  this  type  cf  approximation.  Geometric  relationships 
(sweepback  angles)  and  aspect  ratio  were  preserved  in  the  scaling  up  of 
the  vertical  tail. 

In  addition,  the  cross-sectional  areas  (in  the  body  namelist)  were 
increased  to  include  the  effects  of  the  booms.  Because  of  the  blanking 
effect  of  the  TVU ,  only  the  area  of  one  boom  was  added  to  the  fuselage. 
When  the  TVU  was  absent  from  the  RPV  configuration,  the  areas  of  both 
booms  were  added  to  the  fuselage.  All  other  data  utilised  in  the 
representations  was  taken  from  Table  I. 

Two  parameters  were  chosen  for  investigation:  xcg  location  and 
TVU  attachment.  The  xcg  location  assumed  three  values  (10,  25,  and 
40  percent  mac).  Both  the  TVU  attached  and  unattached  cases  were 
examined  for  each  xcg  location.  Thus,  a  total  of  six  cases  were 
investigated: 


9 


Case  1  -  TVU  attached,  xcg  ~  1 0?S  mac 

Case  2  -  TVU  attached,  xcg  =  25«  mac 

Case  3  -  TVU  attached,  xcg  =  4 DC  mac 

Case  4  -  TVU  unattached,  xcg  =  10%  mac 

Case  5  -  TVU  unattached,  xcg  =  25%  mac 

Case  6  -  TVU  unattached,  xcg  =  40%  mac 

Hence,  the  effect  of  the  TVU  and  the  xcg  location  on  the  RPV  character¬ 
istics  could  be  determined. 

With  the  input  data  finalized,  the  RPV  representations  were  com¬ 
plete.  Figure  2  contains  the  input  data  cards  for  Digital  Datcom. 

The  four  underlined  statements  had  to  be  changed  with  each  case.  A 
visual  interpretation  of  the  RPV  with  the  TVU  attached  and  unattached 
can  be  found  in  Figs  3  and  4,  respectively. 

Data  Reduction 

Since  the  entire  flight  regime  specified  16  Mach  numbers  and  20 
angles  of  attack  per  Mach  number.  Digital  Dotcom  generated  320  sets  of 
stability  derivatives  for  each  case  (i.e.,  a  data  bank).  Unfortunately , 
the  aircraft  equations  of  motion  could  accept  only  one  set  of  stability 
derivatives  at  a  time.  Thus,  a  nominal  flight  condition  was  needed. 
Three  techniques  were  examined  in  the  selection  of  this  flight  con¬ 
dition:  lift  averaging,  best  endurance,  and  best  range. 

The  lift  averaging  technique  assumed  straight  and  level  flight. 

Thus , 

W  =  l 

-  \  pV^SCl 
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3.  Comparison  of  the  Actual  RPV  and  Its  Di 
Da  teem  Representation  with  the  TVU  Attached 


Figure  4.  Comparison  of  the  Actual  RPV  and  Its  Digital 
Datcom  Representation  with  the  TVU  Unattached 
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Substituting  V2  =  M^c2  and  solving  for  Cl  yields 


CL 


_ 2  _W  __ 

SPc2M2 


Assuming  the  maximum  weight  (250  lb)  and  substituting  the  appropriate 
values  gives 


or 


cl 


5.45  x  IQ’3 
M~2 


•  30  _x  JOli 
CLV2 


(1) 

(2) 


From  Eq  (1),  Cl  can  be  calculated  over  the  range  of  reasonable  Mach 
numbers  (Table  II).  A  linear  average  of  lift  coefficient  was 
determined  as  Cl  =  0.59  which  represented  the  average  amount  of  lift 
required  over  the  entire  flight  regime.  Substituting  this  value  into 
Eq  (2)  produced  an  associated  Mach  number  of  Mach  0.10.  Digital  Datcoin 
specifies  its  flight  condition  in  terms  of  angle  of  attack  rather  than 
the  lift  coefficient.  Therefore,  the  data  bank  v/as  interpolated  (using 
Mach  0.10  and  Cl  =  0.59)  to  produce  an  associated  angle  of  attack  of 
a  =  2  deg. 

As  alternate  methods  for  selecting  a  nominal  flight  condition, 
best  endurance  and  best  range  Flight  profiles  were  considered  (Ref  7: 
Sec  III,  33,  37).  After  considering  that  the  RPV  would  be  flown  oy 
sight,  it  became  apparent  that  neither  of  these  two  flight  profiles 
could  be  maintained.  Tnerefore,  the  results  of  the  lift  avoraginu 
technique  (Mach  0.10,  m  =  2  deg)  were  used  as  the  nominal  flight 
condi tion. 
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Finally,  the  trim  condition  needed  to  be  verified  to  insure  that 
this  was  an  appropriate  choice  of  flight  condition.  During  straight  and 
level  flight,  the  pitching  moment  equation  becomes 

Cm  =  Cing  +  Cn1a  a  +  Cm5e  6e 

In  a  trimmed  aircraft,  Cm  =  0.  Using  this  value  and  solving  for  ^etrl-in 
produces 

•  -  "Cfn0  ”  Cm,,  a 

<5ctrnn  "  - a — 

Cr,i6e 


Substituting  the  Case  2  values 


a e trim 


-0.06  +  (0.023) (2 .0) 
'  -0.034 


=  0.33  deg 


Thus,  it  was  decided  that  the  RPV  was  in  a  trimmed  flight  condition. 
Table  III  contains  a  complete  summary  of  the  stability  derivatives  at 
the  nominal  flight  condition. 


Ha_nd  Calculations  of  the.  Remaining  Terms 

Three  types  of  parameters  were  still  needed  to  complete  the  model 
of  the  RPV:  the  remaining  stability  derivatives  (e.g.  CXll),  the  control 
derivatives  (e.g.  Cn^  ) ,  and  the  physical  quantities  (e.g.  Iy) .  The 
complete  derivation  of  these  parameters  is  given  in  Appendix  A.  In 
addition,  these  parameters  are  summarized  in  Table  III. 


Summary 


Table  III  contains  a  summary  of  the  mathematical  models  used  in 
this  study.  Five  approximations  were  included  in  the  generation  of 
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nd  control  derivatives  have  units  of  per  rad. 


these  models: 


1.  Circular  cross-sections  were  used  for  the  fuselage. 

2.  A  straight  tapered  wing  was  substituted  for  the  cranked  wing. 

3.  One  vertical  tail  of  increased  area  (75%  greater)  was  sub¬ 
stituted  for  the  twin  vertical  tails. 

4.  The  effects  of  the  booms  were  incorporated  into  the  fuselage 
representation. 

5.  Average  stability  derivatives  were  computed  using  the  root  wing 
section  and  tip  wing  section  cases. 
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III.  Analysis  of  the  RPV 


Before  any  reasonable  FCS  synthesis  could  be  undertaken,  flight 
specifications  had  to  be  developed  and  the  unaugmonted  system 
analyzed.  Thus,  the  purpose  of  this  chapter  is  to  present  specific 
RPV  handling  qualities,  followed  by  the  above  mentioned  analysis. 

This  analysis  was  divided  into  four  parts:  derivation  of  the  aircraft 
equations  of  motion,  examination  of  the  modal  characteristics,  genera¬ 
tion  of  transfer  functions  and  time  responses,  and  development  of  con¬ 
clusions. 

Handl i ng  Qua!  i  ties 

A  literature  review  resulted  in  only  one  publication  (Ref  14)  which 
specifically  addressed  RPV  handling  qualities;  unfortunately,  sections 
of  this  publication  are  incomplete.  However,  since  all  publications 
reviewed  relied  heavily  upon  the  criteria  set  forth  in  various  Military 
Specifications,  it  was  decided  that  a  comparison  study  of  the  following 
four  specifications  would  provide  the  most  appropriate  source  of  handling 
qualities: 

1.  RPV  Flying  Qualities  Design  Criteria,  AFFDL -TR-76-125  (Ref  14). 

2.  Flying  Qualities  of  Piloted  Aircraft,  MlL-F-o78bl3  (Ref  4). 

3.  Flying  Qualities  of  Piloted  V/STOL  Aircraft,  MIL-F-S33C0  (Ref  5). 

4.  Flight  Control  Systems  -  Design,  Installation,  and  Test  of 
Piloted  Vehicles,  MIL-F-9490D  (Ref  17). 

When  present,  the  specifications  from  AFrDL-TR  -76-1 25  were  used 
since  this  report  specifically  addressed  RPVs.  The  handling  qualities 
in  the  Military  Speci fications  were  used  to  complement  those  areas 
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which  AFFDL-TR-76-1 25  did  not  cover.  Table  IV  contains  the  results 
of  this  comparison  study.  It  consists  primarily  of  eigenvalue  specifica¬ 
tions  which  establish  the  allowable  range  of  values  for  the  roots  of 
the  characteristic  equation. 

Several  restrictions  in  determining  these  handling  qualities  were 
imposed  from  the  onset: 

--  The  general  specifications,  such  as  "Each  rr.ode  will  have 
sufficient  damping",  are  not  repeated  in  Table  IV. 

--  Since  the  aircraft  in  this  study  was  a  RPV,  all  stick  force, 
pedal  force,  and  control  force  specifications  were  omitted.  It  was 
assumed  that  the  servo  would  generate  sufficient  power  to  deflect  all 
control  surfaces. 

--  Since  this  RPV  was  designated  as  a  research  vehicle,  it  was 
assumed  that  it  would  be  flown  on  a  "good"  day.  Therefore,  atmospheric 
disturbance  specifications  were  deleted. 

--  This  RPV  will  be  flown  by  sight;  thus,  display  specifications 
were  omitted. 

--  It  was  assumed  that  the  data  link  'was  sufficiently  reliable 
and  capable  of  transmitting  all  necessary  information. 

--  Redundancy  and  failure  rates  were  not  considered. 

Equations  of  Motion 

The  aircraft  equations  of  motion  are  a  set  of  nonlinear,  coupled 
force  and  moment  equations  which  completely  describe  the  motion  of  an 
aircraft  in  inertial  space.  By  assuming  perturbations  about  straight 
and  level  flight,  those  equations  can  be  decoupled  into  a  longitudinal 
set  and  a  lateral -directional  set.  Linearizing  these  equations  about 
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Table  IV.  Handling  Qualities 


CTassi fication:  Class  III  -  medium  maneuverabi 1 i ty 


Flight  Phase:  Category  A  -  rapid  maneuvering 


Level:  Level  1  -  normal  operation 


Longitudinal  Case: 


tsp 

unSp 

;P 

Residual  Oscillation 
Dynamic  Oscillation 


-  0.35  <  csp  <1.30 

-  See  Fig  5 

-  tp  >  0.04 

-  <  ■  0.60  deg 

-  A  short  term  dynamic  oscillation 
produced  by  a  pulse  elevator 
shall  not  diverge  faster  than  a 
time  to  double  of  15  sec. 


!  Lateral -Directional  Case: 

,:DR 

I  w'iDR 

:  Dutch  Roll  Time  Parameter 

i  „ 

Spiral  Mode 

l 

J 

Roll  Effectiveness 

t 

Roll -Spiral  Coupling 


-  'OR  >0.19 

-  '-npR  >  0.40  rad/sec 

-  '  DR'°n[)R  >  0.35  +  Ar,DR,.,nDR 

where  n0R  =  0-014 

f  i'DR2  i:'/l?'DR  '  20 J 

-  Tr  <  1.4  sec 

-  For  a  steady  state  disturbance 
in  :  of  20  deg,  the  time  to 
double  must  be  greater  than 

20  sec. 

-  >  30  deg  y  in  1.5  sec 

-  None 


an  arbitrary  equilibrium  condition  permitted  the  perturbation  equations 
of  motion  to  be  formed.  These  equations  were  written  in  state  equation 
form  so  that  individual  transfer  functions  and  modal  characteristics 
could  be  determined.  The  complete  derivation  of  the  perturbation  state 
equations  is  found  in  Appendix  3. 

Modal  Analysis 

The  models  of  Chapter  II  wore  used  in  the  perturbation  equations 
of  motion,  resulting  in  twelve  sets  of  state  equations,  six  for  each 
of  the  longitudinal  and  lateral -directional  cases.  Recall  that  a  :,1  itn 
equation  is  of  the  form 

x  =  A  x  +  3  u 

where  the  A  matrix  determines  the  characteristic  equation  of  the  system. 
Each  A  matrix  was  used  in  a  <;  craw  iter  program  (Ref  16)  which  determined 
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the  characteristic  equation,  and  hence,  the  nodal  characteristics  of’ 
each  system.  Table  V  summarizes  the  character! stic  equations  while 
Table  VI  summarizes  the  modal  characteristics. 

A  comparison  of  Table  VI  with  Table  IV  indicated  that: 

--  r,Sp>  <  n^R ■»  TR  satisfied  their  respective  handling  qualities 
by  a  wide  margin. 

--  Cp  satisfied  its  specification,  but  only  by  a  small  margin. 

--  (,;nsp  did  not  satisfy  its  handling  quality  in  cases  4  and  5,  and 
was  at  its  upper  limit  in  the  other  cases. 

--  cqr  and  T$  did  not  satisfy  their  respective  sped  fications . 

--  Since  all  the  longitudinal  roots  were  stable,  no  dynamic 
oscillations  exceeded  a  time  to  double  of  15  sec. 

--  No  roll -spiral  coupling  existed  since  each  mode  had  a  real 

root. 

--  Frequently,  residual  oscillations  are  th*  result  of  nonlinear¬ 
ities  in  the  control  system  (e.g.  dead  zone).  Hence,  no  simple  way 
to  evaluate  these  effects  existed  within  the  scope  of  this  study. 

Two  handling  qualities  required  more  investigation  than  the  above 
comparison:  dutch  roll  time  parameter  ( typ  ngq)  and  roll  effectiveness. 
The  dutch  roll  time  parameter  handling  quality  is  given  in  terms  of  the 
phi  to  beta  ratio  (  |.?/B]or)  .  This  ratio  may  be  found  by  determining 
the  magnitude  of  the  phi  and  beta  envelopes  at  an  instant  of  time  from 
the  rudder  pulse  time  responses  (Ref  4:673).  Utilizing  this  procedure 
and  the  time  responses  from  Appendix  C,  the  dutch  roll  time  parameter 
specifications  were  determined  using  the  average  l^/floR  with  the 
following  results: 
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Table  V.  Characteristic  Equations  of  the  Basic  RPV 


Longitudinal  Case: 

Case 

Characteristic  Equation 

1 

(s2  +■  6.34s  +  39 . 35)  ( s2  +  0.040s  t-  0.16)  =  0 

2 

(s2  +  7.38s  +  37.01 )(s2  +  0.042s  +  0.15)  =  0  j 

1 

3 

(s2  +  8.20s  +  24.16) (s2  +  0.052s  +  0.11)  =  0  1 

4 

(s2  +  9.64s  +  65 . 7  2 ) ( s 2  +  0.040s  +  0.18)  =  0 

5 

(s2  +  12.73s  +  31 .49) (s2  +  0.046s  +  0.16)  =  0 

6 

(s  +  11 .27) (s  +  4.68)(s2  +  0.074s  +  0.093)  =  0 

| 

j  Lateral -Directional  Case: 

1 

Case 

Characteristic  Equation 

1 

(s  +  11 . 50) (s  -  0.042) (s2  +  1.49s  +  19.03)  =  0 

2 

(s  +  11 . 1 8 ) ( s  -  0.042) (s2  +  1.42s  +  17.76)  =  0 

3 

(s  +  10.88) (s  -  0.042) (s2  +  1.34s  +  16.50)  =  0 

4 

(s  +  1 1 .30) (s  -  0.044) (s2  +  1.23s  +  17.85)  =  0 

5 

(s  +  10.99) (s  -  0.043) (s2  +  1.17s  +  16.63)  =  0 

6 

(s  +  10.71 )(s  -  0.044)(s2  +  1.08s  +  15.41)  =  0 
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Table  VI.  Modal  Character! sties  of  the  Basic  RPV 


JV-'lDR 

1.63 
1 .87 


t 

i 

i 


( 


Case 

2 


lljr.dl  ing  Qua]  i  ty 
■'-DR'-npR  51  0-^8 
^DR'-  npR  >  0-78 


Comparing  these  values  with  Table  VI  showed  that  the  dutch  roll  time 
parameter  specification  was  satisfied. 

Likewise,  the  roll  effectiveness  of  the  RPV  needed  to  be  determined. 
With  maximum  deflection  of  the  ailerons  (15  deg),  the  RPV  produced  62 
deg  and  64  deg  of  bank  in  1.5  sec  for  cases  2  and  5,  respectively  (Fig  6). 
Hence,  the  RPV  had  over  twice  the  required  roll  effectiveness. 

Although  not  specifically  mentioned  in  the  handling  qualities,  one 
other  characteristic  needed  attention,  the  ratio  of  the  oscillatory  com¬ 
ponent  to  the  average  component  of  bank  angle  ( ;osc/;ave) .  This  ratio 
is  determined  by  analyzing  the  peaks  of  a  bank  angle  time  response  due 
to  an  impulse  aileron  command.  Since  r,DR  <  0.20,  the  following  empirical 
formula  could  be  used: 


Jose  =  *1  _+_  :  3  "  z'-2 
Vve  '.b"i  +  W-j  T_2g" 


where  fq  ,  f 2 »  and  *3  corresponded  to  the  first  second,  and  third  peaks 
(Ref  xxii).  Hence,  the  impulse  ciine  responses  (Fig  7)  were  used  in 
Eq  (3)  and  produced 

Case  fosc/fave 

2  -0.0081 

5  -0.0108 


Thus,  the  RPV  was  not  very  oscillatory. 
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Ira  ns  for  Functions  and  Time  Responses 

before  transfer  functions  and  time  responses  could  be  generated, 
the  actuator  dynamics  must  be  included  in  the  equations  of  motion. 
Modified  Futaba  FPS-14  servos  were  used  for  the  RPV's  actuators.  Test 
of  these  servos  produced  two  different  representations: 


TF  _  (15^7)2 _  _ 

1  hservo  s2  +  ("2T(VdT(T5 . 7)s  +  TT5.7~}2 

(14.45)2 

TFservo  -  STT“(T)‘(78)0 4 . 4 5 ) s  +  (1 4.45) 2 


(4) 

(5) 


Equation  (4)  corresponded  to  a  5.5  in. -lb.  load  and  5.0  volt  input. 
Equation  (5)  corresponded  to  a  11.8  in. -lb.  load  and  10.0  volt  input. 
Since  the  RPV  would  have  a  peak  input  of  22  volts,  some  type  of  extra¬ 
polation  of  Eqs  (4)  and  (5)  was  needed.  It  was  decided  that  a  50  percent 
reduction  in  break  frequency  in  Eq  (5)  would  adequately  compensate  for 
the  larger  input  and  whatever  slippage  that  occurred  due  to  the  linkage. 
Hence,  the  following  representation  was  utilized  for  the  RPV  actuators: 


TF 


(/.22)2_  _ 

servo  -  sZ"+  {271787(7 .227s  +T7722T2 


52.20 

s?  +  Tl.5'6  s  t-  52.20 


(6) 


The  corresponding  differential  equation  is 


6*  -  -52.204  -  11.56?  +  S2.204cm(j  (7) 

The  equations  of  motion  were  augmented  with  Eq  (7)  resulting  in  the 
open  loop  state  equations  given  in  Fig  8.  From  these  equations,  the 
transfer  functions  for  the  basic  RPV  were  generated  (Ref  16).  Only  Case 
2  is  presented  in  this  report  (Table  VII)  since  the  FCS  design  will  be 
based  on  this  case. 
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-0.05S  0.33  -0.29 


Table  VII.  Transfer  functions  tor  the  basic  P.fV,  Case  2 
Longitudinal  Case: 

=  { s 2  +  7.33  s  +  37.01  )(s2  +  0.042  5  <  0.15)(s:-  +  11. 66  s  ->  32. 20) 

=  -2.34  (s  +  4. 76) (s  -  163.03) 

=  -19.04  (s  +  11 4. 58) ($2  +  0.57  s  +  0.28) 

-  -2219.32  (s  +  0.089) (s  +  2.32) 

=  -2219.32  s ( s  +  0 . 089 ) ( s  +  2.82) 

=  52.20  ( s2  +  7.38  s  +  37.01)(s2  +  0.042  s  +  0.15) 

=  52.20  s(s2  t  7.33  s  +  37.01)(s2  +  0.042  s  +  0.15) 

Lateral -Directional  Case: 

v  =  (s  +  11 . 18) (s  -  0.042) (s2  +  1.42  s  +  17.76) (s2  +  11.56  s  +  52. 20)2 

N^cmd  =  -1-24  +  °-26)(s  '  36-7°) 

=  1.23  ( s  -  2 .27) (s  +  2.54)  (s  -  26.33) 

ucmd 

Njcmd  =  43.33  (s2  t  1.43  s  +  15.19) 

M'i^cmd  =  43'83  s  +  1  -43  s  +  15-19) 

fTa  -  52.20  (s  +  1 1 .18)  (s  -  0 . 042 )  ( s2  +  1.42  s  +  !7.76)(s2  *  II.  Co  s  *  :) 

acmd 

N.1;  d  -  52.20  s(s  +  1 1 .18)  (s  -  0.042)(s2  +  1.42  s  +  17.76)(s7  t  11. 5C  s  +  32.20) 

?i*r  -  o 

'  c1cmd 

-  0 

acmd 

.V?cnid  =  °-51  (s  -  0.073)  (s  +  11  .91 )  (s  +  79.23) 

Nr  -•  -40.04  (s  +  11.65)(s2  +  0.24  s  +  0.49) 

rcmd 

N .•  rcmcj  =  13.67  (s  +  4.32)(s  -  14.21) 

NPrcmd  =  13-67  s(s  +  4 . 32) (s  -  14.21) 

N  rcmd  =  0 

"Ted  ' 0 

N.^md  =  32.20  (s  +  11  . 1 8 ) ( s  -  0.042)(s2  +  1.42  s  t  17.76)(s2  +  ]]  .  or,  t  ',7.20) 

-=  52.20  s(s  +  11  . 18) (s  -  0.042)  (s2  «•  1.42  s  +  17.76)  (s2  >  11.56  s  *  57. ;n) 


N 


N 


ecmd 

a 

f  ecmd 
9 

!'ecmd 

NKmd 

r.e 

ecmd 

N  "e 
v  ecmd 


fl  r 
rcind 


Typically,  a  pilot's  command  consists  of  an  input  and,  a  short 
time  later,  removal  of  that  input.  With  this  in  mind,  a  1  see  pulse 
was  selected  as  the  input  for  all  time  responses.  The  magnitude  of 
this  pulse  was  -5  deg  for  the  elevator  command  and  i5  deg  for  the 
aileron  and  rudder  commands .  A  set  of  time  responses,  for  case  2,  were 
generated  (Ref  16)  using  these  inputs  and  are  contained  in  Appendix  C. 

Even  though  only  one  input  was  used,  several  general  observations 
were  made. 

--  All  the  time  responses  had  the  shape  that  would  be  expected 
far  this  input  (Refs  2:47,  120-130  and  12:306,  668). 

The  peak  values  were  high.  This  stormed  from  the  slow  airspeed 
and  large  control  surface  effectiveness. 

--  The  phugoid  nniural  frequency  was  much  higher  than  most 
aircraft. 

--  The  dutch  rail  mode  was  very  oscillatory. 

--  The  spiral  "Oiie  had  a  dominant  effect  that  could  not  he  ignored. 
Cone  ]  nr.  ions 

Several  conclusions  wore  drawn  fro;:  the  modal  analysis  and  time 
res nr-e .  Th  -so  conclusions  define  the  design  objectives  for  Chapter 
IV.  In  the  1'V-ji  turlinal  rase,  nhugn  j.!  natural  frequency  needed  to  be 
decreased  r i nee  oscillations  of  tin's  mennilur/'  and  frequency  wore 
undesirable.  In  the  1 atoral -d i rec t ion  al  case,  the  effect  of  dutch 
roll  oscillation  needed  to  he  reduced  and  the  spiral  mode  needed  to 
be  stabilized. 

Tiie  effect  of  the  TVU  on  the  RPV  was  minimal.  All  natural 
frequencies  vere  slightly  reduced  as  was  the  short  period  damping 
ratio.  ihc  >’  ,er  characteristics  remained  essentially  unehinged. 


33 


A:,  the  xcg  location  move d  further  aft,  the  following  characteristics 
were  increased:  .sp,  p,  and  .  n^.  Also,  Urn  was  decreased.  All 
other  characteristics  regained  essentially  unchanged. 
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IV.  Res  inn  and  Eva! nation 

All  design  processes  consist  of  the  same  procedure:  synthesize, 
evaluate,  and  resynthesize.  The  purpose  of  this  chapter  is  to  complete 
this  design  process  for  a  PCS  to  correct  the  deficiencies  found  in 
Chapter  III.  The  synthesis  phase  consisted  of  a  longitudinal  design 
and  a  lateral -directional  design.  The  evaluation  phase  consisted  of 
three  sensitivity  studies,  a  final  handling  qualities  check,  the  gen¬ 
eration  of  transfer  functions  and  time  responses,  and  a  nonlinear 
simulation. 

L  o  ng  i  _t 1  id  trial  Dos  i  cj  n 

Even  though  the  longitudinal  case  satisfied  all  its  handling  qual¬ 
ities,  the  phugoid  mode  was  quite  oscillatory.  Thus,  it  was  desirable 
to  reduce  this  oscillatory  response  oiteer  by  increasing  p  or  de¬ 
creasing  !,np.  Because  of  the  constraint  of  only  rates  for  acceptable 
feedbacks,  the  choices  of  feedbacks  for  augmentation  wore  pi tch  rate 
ana  pir.cn  *ui l.e  with  compensation. 

It  -  'ui  red  only  a  brief  glance  at  the  root  locus  of  trie  pitch 
rare  fee-abac k  sys t.-ra  (Pig  9)  to  realize  the t  it  would  not  provide 
appropriate  augi .  z : ,  La  t i o  n .  This  was  due  to  the  zero  that  was  very  close 
to  the  phugoid  poles  whim  tr.'v;  lato.i  into  a  vary  large  gain  for  a  small 
increase  in  camping  In  its- If,  the  large  gain  was  not  a  problem; 

rather,  it  was  its  effect  on  the  other  roots  that  presented  the  diffi¬ 
culty.  A  large  gain  would  cause  the  short  period  roof,  to  become  less 
stable  and  its  handling  quality  to  become  unsatisfied.  In  addition, 
the  actuator  dynamics  would  change  gr  eatly  which  •/;•••  s  ne-iesi r  ihle . 
Therefore',  tne  conclusion  was  t'-at  the  addition  of  mi  me nsub’Vn  would 
be  advantageous. 
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The  desire  to  keep  the  compensation  si;  pic  and  to  eliminate  the 
effect  of  the  zero  at  the  origin  dictated  the  choice  of  a  pure  integrator 
as  the  feedback  compensator .  The  root  locus  (Fig  10)  revealed  that  the 
zero  at  the  origin  was  effectively  removed  and  that  the  short  period 
root  b  ;camo  unstable  faster  than  without  compensation.  The  short 
period  handling  qualities  restricted  the  amount  of  feedback  gain  to  less 
than  5.14.  A  feedback  gain  of  0.050  produced  the  best  compromise  between 
maximum  augmentation  and  minimum  effect  on  the  other  roots.  Thus,  the 
control  law  was 

;;f?cmd  (s)  '•  Fin  (s)  f  A.-PJ?P.  q  ( s )  (8) 

An  examination  of  the  phngoid  characteristics  showed  that  the  damping 
ratio  increased  from  -,p  -  0.055  to  op  =  C.25  while  the  natural  frequency 
remained  essentially  constant.  Hence,  the  oscillatory  phugoid  mode  was 
reduced  substantially. 

Lateral  directional  Design 

The  analysis  of  Chapter  III  suggested  that  the  dutch  roll  damping 
ratio  should  be  increased  and  the  spiral  mode  should  be  stabilized.  A 
standard  yaw  damper,  consisting  of  yaw  rate  feedback  to  the  rudder 
(passing  through  a  washout  circuit),  was  used  to  augment  the  dutch  roll 
mode.  The  main  difficulty  in  this  design  was  selecting  the  washout  pole 
location.  A  large  time  constant  for  this  pole  indicates  a  large  in¬ 
crease  in  both  damning  ratio  and  uncoordinating  effect  (Ref  ]?:(•  1 )  • 

A  washout  pole  location  of  s  =  -1.75  (T  -  0.57  sec)  produced  the  lest 
compromise  to  the  trade-off  described  above.  A  feedback  gain  of  0.033 
was  selected  because  it  corresponded  to  I  he  maximum  dutch  roll  damning 
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ratio  (see  Fig  11).  Thus,  the  control  lew  was 

' rCind  ~  ^in  (s)  +  s  +  ]  '75  r  (9) 

An  examination  of  the  dutch  roll  characteristics  revealed  11. at  the 
damping  ratio  increased  from  "dr  -  0.17  to  -  0.29  while  the  natural 
frequency  increased  from  -mtqr  =  A. 21  rad/sec  to  .. nrjp,  =  5.12  rad/sec. 

To  simplify  the  system,  the  v:ashout  circuit  was  removed  from 
E  ,9).  With  this  control  law,  a  feedback  gain  of  0.032  produced  the 
n  .‘u:n  dutch  roll  damping  ratio  (see  Fig  12).  The  dutch  roll  char- 

a>. , eristics  revealed  that  the  damping  ratio  increased  to  =  0.21. 
Unfortunately,  this  increase  was  not  enough  to  satisfy  the  handling 
qualities  by  an  adequate  margin.  Thus,  the  washout  circuit  was  rein¬ 
serted  into  Eg  (9 ) . 

The  spiral  mode  was  most  dominant  in  the  lateral  tine  responses. 
Thus,  either  integral  of  roll  rate  or  yaw  rate  feedback  to  the  ailerons 
would  be  the  most  effective  means  of  stabilizing  the  spiral  mode.  For 
the  integral  of  roll  rate  feeback  system  (fig  13),  a  feedback  gain  of 
0.012  brought  the  spiral  root  just  inside  the  left  hand  plane  resulting 
in  a  control  law 

■‘•cmd  (0  =  Ain  (s)  +  p  (s)  (10) 

Implementation  of  Fq  (10)  had  no  effect  on  ar.y  lateral-directional 
characteristic  except  the  spiral  root  which  moved  to  s  =  -0.0003  from 
s  -  0.041 . 

For  the  yaw  rate  feedback  system  (Fig  14),  a  feedback  gain  of 
0.044  established  0  stable  spiral  node, 
system  was 

3R 


The  control  law  for  this 


Figure  11.  Root  Locus  for  the  Yaw  Damper  Including  the  Washout  Circuit 


a 


Figure  12.  Root  Locus  for  the  Yaw  Damper  Excluding  the  Washout  Circuit 
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Figure  14.  Root  Locus  for  Yaw  Rate  Feedback 
to  the  Ailerons 
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^acmd  (s)  ~  ^in  (s)  ~  0.044  r  (s) 


(11) 


Implementation  of  Eq  (11)  also  had  no  effect  on  any  lateral -directional 
characteristic  other  than  moving  the  spiral  root  to  s  =  -0.0001  from 
s  =  0.041. 

Thus,  both  Eqs  (10)  and  (11)  accomplished  the  desired  goal  of 
stabilizing  the  spiral  mode.  Equation  (11)  was  selected  as  the  more 
desirable  control  law  because  it  did  not  require  any  compensation  or 
any  additional  instrumentation,  such  as  a  roll  rate  gyro. 

Sensi ti vi t.y  Studies  and  Redesign 

Three  types  of  sensitivity  studies  were  used  to  evaluate  the 
PCS:  a  flight  condition  study,  a  parameter  study,  and  a  physical 
quantity  study.  Each  study  based  its  results  on  a  modal  analysis. 

The  flight  condition  study  varied  the  airspeed  throughout  the  entire 
speed  regime.  The  parameter  study  varied  each  parameter  in  the 
equations  of  motion  (except  Uo»  q,  S,  and  c)  by  :*  25  percent.  Lastly, 
the  physical  quantity  study  varied  the  vertical  tail  area. 

The  modal  analysis  from  the  flight  condition  study  showed  that 
three  handling  qualities  were  not  satisfied  over  the  entire  airspeed 
regime.  First,  the  short  period  natural  frequency  handling  quality 
became  unsatisfied  at  the  higher  airspeeds  (M  i  0.12).  Feedback  gain 
adjustment  had  no  effect  on  increasing  the  range  of  airspeeds  that 
satisfied  this  handling  quality.  Thus,  a  restriction  was  imposed 
stipulating  that  this  handling  quality  would  not  be  satisfied  at  the 
higher  airspeeds.  This  restriction  was  acceptable  when  considering 
that  the  higher  airspeeds  were  clearly  in  the  upper  limits  of  the 
RPV's  capabilities.  Secondly,  this  study  revealed  that  the  spiral 
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inode  was  too  unstable  at  low  airspeeds.  Increasing  the  feedback  gain 
in  Eq  (11)  to  0.09  stabilized  this  mode  throughout  the  entire  airspeed 
regime.  Lastly,  this  study  found  that  the  dutch  roll  damping  ratio 
did  not  satisfy  its  specification  at  the  higher  airspeeds  (M  -  0.15). 

A  feedback  gain  of  0.03  maximized  this  dampino  ratio  throughout  the  air¬ 
speed  regime.  Unfortunately,  as  in  short  period  natural  frequency,  the 
entire  airspeed  regime  could  not  be  satisfied  and  the  higher  airspeeds 
restriction  was  imposed.  Again,  this  was  acceptable  when  considering 
the  RPV's  capabilities. 

Hence,  the  control  laws  (Kgs  (3),  (9),  and  (11))  were  refined  to 

Secmd  (s )  r  Ejn  (s)  +  q  (s)  (12) 

'iacmd  (si  "  A (n  (s)  -  0.090  r  (s)  (13) 

firemd  (s)  =  Rin  ( s )  +  !?-'  °/y  75  r  (s)  (14) 

Plots  of  the  pole  locations  due  to  a  varying  airspeed,  with  the  above 
control  laws  implemented,  may  be  found  in  Fig  15.  These  plots  show 
pole  locations  for  airspeeds  from  Mach  0.05  tc  Mach  0.15  at  Mach  0.01 
intervals.  The  arrows  indicate  increasing  airspeed. 

The  second  study  entailed  tabulating  the  modes  of  the  RPV  as  each 
parameter  in  the  equations  of  motion  (except  Uq,  q,  S,  and  c)  was 
varied  u25  percent.  Comparing  the  extremum  (Table  VIII)  with  the 
handling  qualities  (Table  IV)  showed  that  all  handling  qualities  were 
satisfied.  Thus,  an  error  in  any  one  parameter  of  less  than  25  percent 
will  not  cause  any  handling  quality  to  be  exceeded.  Of  course,  r.o 
assurance  can  he  made  if  mote  than  one  parameter  changes. 
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Table  VIII.  Extremum  of  the  Parameter  Sensitivity  Study 


Source  of  the  extremum  is  given  in  parentheses.  The  sign 
indicates  whether  the  parameter  was  increased  or  decreased 


1 


The  last  study  was  completed  because  of  the  single  vertical  tail 
approximation.  Three  scale  factors  (1.50,  1.75,  and  2.00)  were 

examined  to  determine  the  effect  of  the  vertical  tail  area  on  the  system 

modes.  As  expected,  the  vertical  tail  area  had  no  effect  on  the  long¬ 
itudinal  case.  Little  effect  was  found  on  the  lateral-directional  case. 
A  plot  of  pole  locations  due  to  varying  vertical  tail  area  is  shown  in 

Fig  16.  The  arrows  indicate  increasing  vertical  tail  area.  Comparing 
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Figure  16.  Pole  Movement  due  to  Variations  in  Vertical  Tail  Area 


Fig  16  and  Table  IV  revealed  that  all  handling  qualities  were  satisfied. 
Thus,  a  substantial  variation  in  vertical  tail  area  did  nut  cause  any 
handling  quality  to  become  unsatisfied. 

Handling  Qua!  i  ties  Check_ 

A  final  check  of  tne  handling  qualities  was  made  with  the 
control  laws  (Eqs  (12)  -  (14)).  A  comparison  of  t ho  augmented  RPV, 
the  basic  RPV,  and  the  handling  qualities  (Table  IX)  revealed  that  all 
characteristics  satisfied  their  respective  handling  qualities. 


45 


Table  IX.  Comparison  of  the  Augmc n ted  RPV, 
Basic  RPV,  and  the  Handling  Qualities 


Augmented 

Basic 

Hand 1 i ng 

RPV 

RPV 

Qual i ti es 

?sp 

0.58 

0.61 

0.35  <  rsp  <1.30 

unSp,  rad/sec 

5.74 

6.08 

Fig  5 

tp 

0.25 

0.055 

Cp  >  0.04 

wnp,  rad/sec 

0.41 

0.39 

★ 

Tr,  sec 

0.087 

0.039 

Tr  <  1.40 

TS,  sec *  ** 

22.91 

-24.04 

Ts  >  0  or  Ts  <  -23.85 

(") 

(16.66) 

(>20.00) 

?DR 

0.31 

0.17 

'dr  >0.19 

^noR,  rad/'sec 

4.60 

4.21 

u dqr  >  0.40 

Dutch  Roll  Time 
Parameter,  sec"1 

1  .44 

0.71 

4DR  WnDR  >  0.35 

*  No  handling  quality  exists. 

**  Negative  sign  indicates  instability.  Time  to  double  (sec)  is  given 
in  parentheses. 
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As  before,  the  dutch  roll  time  parameter  (  g;;  /  c j ec i fi cation  was 

dependent  on  the  phi  to  beta  ratio  (.;/•! DR )  -  ‘’he  time  responses  of 
Appendix  C  revealed  that  the  graphical  method  of  determining  ! : /  ■  i (;r  used 
in  Chapter  III  could  not  be  utilized.  Therefore,  an  analytical  method  was 
used.  The  i ; /■■  i DR  is  the  ratio  of  the  magnitude  of  the  .  component  to 
the  magnitude  of  the  r  component  of  the  dutch  roll  eigenvector.  Thus, 


Since 


I  -0^00  7  5  j  0 .001 7 
I  0 . Out* 4  +  J0.0567 

=  0.13 


2.84 


f  20.00 


the  dutch  roll  time  parameter  specification  became  ' gp  '-rtoR  0.3-3, 
which  was  satisfied  by  the  augmented  RPV. 

In  addition,  the  roll  effectiveness  of  the  juijmcntod  RPV  needed 
to  be  computed.  With  maximum  deflection  of  the  ailerons  (15  deg),  the 
augmented  RPV  produced  62  deg  of  bank  in  1.5  oc  (Fig  17).  Hence,  no 
change  in  roll  effectiveness  was  experienced  due  to  the  FOS. 

Lastly,  since  0.20,  Fq  (3)  could  not  be  utiliz  <1  to  determine 
•'osc/°ave-  Rather 


‘ use  _  !1  -  :2 

;JVe  i  r  l 


should  be  used 
was  used  in  Eg 


(Ref  4:  xxii).  Thus ,  the  impulse 


(15)  and  produced 


■  o-.r 
•  ave 


0.044 


time  rt  spouse 


(I  ’g  13) 
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Nonce,  tie  augmented  RPV  was  not  very  oscillatory. 


Funct ions  and  Time  Responses 

Before  transfer  functions  and  time  responses  could  be  determined, 
the  closed-loop  state  equation  must  be  derived.  Recall  the  form  of 
the  open-loop  state  equation 

x  =  A  x  +  8  u  (16) 

where 

u  =  K  x  +  IN  (17) 

Substituting  Eq  (17)  into  Eo  (16)  and  simplifying 

x  =  Ax  +  B(Kx  +  IN.) 

=  (A  +  BK)x  +  B  IN  (18) 


Substituting  the  control  lows  (Cqs  (IB)  -  (la)  into  Lq  (18)  provided 
the  closed-loop  state  equation  given  in  Fig  If*.  Notice  that  a  new  state 
variable  had  to  he  added  to  the  lateral -directional  state  equation  due 
to  the  washout  circuit  in  the  vaw  damper.  The  new  state  was  chosen 
to  be  the  yaw  damper  feedback  B  and  defined  as 


cl)  -  0-030  s  /  > 
s  V  1  .75  (  ' 


Using  the-  closed-loop  state  equation  (Fig  19),  transfer  functions 
and  time  responses  of  tne  augmented  RPV  were  gt  r.era  ted  (Ref  10).  Table 
X  contains  the  resultant  Iran  ;f.rr  functions  while  the  time  responses 
are  in  Appendix  C.  Comparing  thece  time  responses  with  the  basic  RPV 
time  responses  verified  t>  at  the  d.-sign  ebjoctives  had  bn  n  achieved. 


Tabu?  X.  Transfer  functions  for  the  Auucon tod  P.PV 
Lornii tu J i nal  Case: 


1  V 

=  (sT* 

0.21  ? 

>  +  0. 1 7) ( s2  +  6 . 62  s 

+  32.38) (s2  +  12 

=  -2.34 

(s  +  * 

1 .76) (s  -  163.08) 

Mn 

NEin 

=  -19.04 

-  (s  + 

114.58) (s2  +  0.57  s 

+  0.28) 

i  o 

; 

-  -2219. 

32  (s 

+  0 . 089 ) ( s  +  2.32) 

1  N9. 

tin 

=  -2219. 

32  s(« 

;  +  0.089) (s  +  2.82) 

Nr? 

=  52.20 

(s2  + 

7.38  s  +  37.01 )(s2  + 

0.G42  s  +  0.15) 

;  On 

!  N;:e 

-  52.20 

s(s2  -1 

-  7.33  s  +  37.01 ) (s2 

+  0.042  s  +  0.15) 

i  Cin 

Lateral -Directional  Case: 


-  Is  +  11. 50) ( s  +  0.0*14)  (s  +  2.32)(s2  +  2.89  s  +  21.12) 

(s2  +  11.56  s  +  52 . 2  0)4; 2  +  9.11  s  +  31.49) 

---  -64 . 78 (s  +  O.ZT  )(s  +  2.4';) (s  -  3S.72)(s2  +  10.95  s  +  46.42) 

=  G4 . 09 ( s  +  1 .75) (s  -  Z . 27) (s  +  2 .54) (s  -  26.4^1 
(sZ  +  11.56  s  +  62.18) 

=  2283. 02(s  +  2. 30) (s2  +  2.96  s  +  20. 39) (s2  +  9.49  r>  +  29.63) 

-  2283.C-2  s ( s  +  2.30)(s2  +  2.95  s  +  20.39) (s2  +  9.49  s  +  29.63) 

-  52 . 20(s  -  0.041 ) (s  +  2.31  )(s  +  11 .10)  (s2  9.32  s  +  30.75) 

(s2  +  2.63  s  +  23.44) 

=  52.20  s(s  -  0.041 )(s  +  2.31 ){$  +  1 1 .10) (s2  +  9.82  s  +  30.75) 
(s2  +  2.63  s  +  23.44) 

=  100.27  (s  -  0.C0C4 ) (s  -  2.27)(s  +  2.55) (s  -  26.32) 


Tabic  X.  (Continued) 


I  Ain  =  100.27  s(s  -  0.0004)(s  -  2.27) (s  +  2.s5)(s  -  25.32) 

I  ilL  -  1.02  ( s  -  0. 0304)  ( s  -  2 . 2  7 )  ( 5  +  2.55)  (g  -  2G.32) 
i  (s-1  r  11.55  s  +  52.20) 

!  flR.  -  26.45  (s  +  0.019)  (s  -r  1 .  75)  (s  f  7 1 . 74 ) ( g  +  79.24) 

;  '  (s2  +  11.63  s  +  51.16) 

HRin  =  -2090.00  (s  +  1  . 75) Cs  •-  11 . 6 5 )  ( s 2  +  0.24  s  +  0.49) 

(s'-  +  11.56  s  +  52.20) 

Nr j  =  713.53  (s  +  1  . 75) ( s  -  10.47) (s  -  1 2.95) (s2  +  4.15  s  +  22.30) 

U P  =  713.53  3 ( s  +  1 . 75) (s  +  10.47)(s  -  12.95) (s2  +  4.15  s  +  22.30) 

Kin 

=  9313.84  (s  +  1.761(5  f  11  .65) (s-  +  0.24  s  +  0.49) 

j  Ki  n 

i  =  9813.84  s(s  +  1.75)(s  +  1 1 .66) (s2  +  0.2-1  s  ^  0.49) 

I  K  i  n 

I  Npr  52,20  (s  f  0 . 04 6 )  ( s  +  1 .75)  (s  +  1 1 .52)  (s2  +  1.37  s  +  16.24) 
j  ln  (s'-  +  11.19  s  +  53.26) 

Np'r  -  52.2  s(s  <-  0.045)(s  +  !.75)(s  +  ?1.52)(s2  +  1.37  s  +  16.24) 
Kin  ($2  +  n.19  s  !-  53.26) 

Nn.  -  -62.70  (s  -  0.0003) (s  +  1 1 .65) (s2  +  0.24  s  (-  0.49) 
i  Aln  (s2  +  11.56  s  +  52.20) 
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Nonl inear  Simulation 

The  last  form  of  evaluati  an  was  a  complete  nonl  inear  simulation 
which  examined  the  effect  of  coupling  on  the  PCS.  Initially,  an 
analog  simulation  was  completed  with  the  results  being  of  low  quality. 
Therefore,  a  digital  simulation  was  accompl i shed  using  the  numerical 
integration  subroutine  ODE  (Ref  13)  (see  Fig  20).  The  equations  of 
motion,  in  their  most  basic  form  (Eqs  (35)  and  (33)),  were  used  as  was 
the  same  control  inputs.  Recall  the  nonlinear  states  (except  U)  have 
the  same  value  as  the  perturbation  states  since  the  equilibrium  con¬ 
ditions  were  chosen  zero.  Appendix  C  contains  the  results  of  this 
nonlinear  simulation.  A  comparison  with  the  other  time  responses 
indicated  that  the  coupling  effect  between  modes  was  quite  small. 

Thus,  the  decoupling  assumption  produced  minimal  error  and  was  accept¬ 
able. 

In  the  formulation  of  the  simulation  algorithm  (Fig  20),  ore 
critical  selection  was  made,  the  increment  time.  Shannon's  sampling 
theorem  dictated  that  ft  d  0.68  sec  for  complete  reproduction  of  the 
response  (Ref  9:  Sec  III,  21).  In  practice,  a  much  faster  increment 
time  than  At  *  0.68  sec  is  desirable.  Thus,  At  -  .02  sec  was  chosen 
for  the  first  2.0  sec  of  the  simulation  and  At  -  0.10  sec  for  the 


remainder  of  the  simulation. 
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V.  A  Practical  /Application 

This  study  would  not  be  complete  without  addressing  an  application 
of  this  RPV.  It  was  hoped  that  by  demonstrating  its  usefulness,  motiva¬ 
tion  for  continued  research  in  this  area  would  be  generated.  The  chosen 
application  was  termed  model  matching,  i.e.,  making  the  RPV  have  similar 
dynamic  characteristics  as  another  aircraft  through  the  design  of  an 
expanded  FCS.  The  aircraft  chosen  for  duplication  was  the  F-15.  The 
first  step  of  this  procedure  was  to  determine  the  F-I5's  dynamic  char¬ 
acteristics.  Then,  three  design  procedures  were  investigated  for 
accomplishing  this  duplication:  Gui llemin-Truxal  design.  Entire 
E i gens  true turo  Assignment,  and  an  extension  of  Observer  theory.  Only 
the  longitudinal  case  was  examined  since  the  lateral -directional  case 
would  be  a  repetition  of  the  same  theory. 

Analysis  of  the  F-15 

The  mathematical  model  for  the  F-15  was  obtained  from  Ref  19  and  is 
contained  in  Table  XI.  As  with  the  RPV,  the  equations  of  motion  were 
augmented  with  the  actuator  dynamics.  The  following  representation 
was  utilized  for  the  F-15  servo  (Ref  19): 

—  20 
'f servo  "  s  VT0 

The  corresponding  differential  equation  was 

6e  -  -20  .e  1  20  S Ocmd  (19) 


The  equations  of  motion  were  augmented  with  Eg  (19)  resulting  in  the 
longitudinal  state  equation  given  in  Fig  21.  Transfer  functions 


lab  1  e  Xl.  F - 1  5  Mdtr.c-ii.d tical  Model 


m  =  10B5  sluqs 
Uq  =  665  ft/sec 
q-  =  333  lb/ ft2 
S  =  608  ft2 
c  =  16  ft 

Iy  -  155746  si ucj-f t2 
Cxu  =  --04 
Cx  •  =  0 

*rt 

Cx  -=  0.24 

ot 

Cxq  =  0 
cx  >e  “  0 
C2u  =  -0.30 

C|_.  =  i.oi 

CLa  -  3.78 
CL^  -  3.12 
c2>5e  =  0 

C|,i(j  =  0 

Cm-  =  -1.20 
Cmq  =  -3.70 
Cm.ce  -  -0.69 


stability  and  control  derivatives  have  the  units  of 
rad . 
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(Table  XII)  and  time  responses  (Figs  25  -  29)  were  generated  from  this 
state  equation  (Ref  16).  Also,  the  longitudinal  characteristics  were 
determined  as 


(~nSp  =  2.90  rad/sec 

GP  -  0.10 

uOp  =  0.063  rad/sec 
Gui 1 1 emi n-Truxal  Design 

The  Guillemin-Truxal  design  procedure  is  basically  a  pole/zero 
placement  technique.  From  a  desired  closed-loop  transfer  function, 
the  appropriate  cascade  or  feedback  compensator  is  mathematically 
calculated  (Ref  6:408).  This  procedure  guarantees  duplication  of  the 
desired  closed-loop  transfer  function  for  single  input-single  output 
systems.  However,  in  this  application,  the  RPV  was  a  single  input 
Usciiid) _|1HJltiple  output  ('u,  «,  o,  q,  -:e)  system;  therefore,  the 
validity  of  this  technique  was  questionable.  However,  it  was  decided 
that  the  duplication  of  one  state  might  prove  valuable.  Of  course, 
the  duplication  of  one  state  would  alter  the  other  three  states  which 
might  provide  a  closer  match  to  all  of  the  F - 1 5  states. 

For  a  cascade  compensator  Gc(s),  the  closed-loop  unity  feedback 
transfer  function  is 

N(s)  _  G r ( s )  G ( s ) 

J(s)  "  1  +  Gc (s)G(sT 

Solving  for  Gc(s) 

r  f  c  )  -  H  (  5  ) 

Gc(')  ■  '['D ( s )  -  n(,)]G(s) 
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Applying  this  equation  to  the  RPV  .and  E-lb  st cites  (Ref  16)  r.jsul  ted  in 
four  cascade  compensators  (Table  XIII),  all  of  which  were  very  compli¬ 
cated.  From  a  hardware  point  of  view,  these  compensators  would  be 
extremely  difficult  to  implement,  if  not  impossible.  Thus,  cascade 
compensators  were  rejected. 

This  same  procedure  was  repeated  for  feedback  compensators.  The 
closed-loop  transfer  function  is 


N(s)  _  G_(s) 
ii'vs)  "  i  +  bcTT^CsT 


Solving  fer  G^s) 


Gc  ( s ) 


1)  (  S  )  G  ( S )  -_Njs] 
hs)G(s) 


Tilis  equation  resulted  in  four  feedback  ccmpensa tors  (Table  XIV)  which 
were  simpler  than  the  cascade  compensators,  but  still  impractical  to 
implement.  Since  botii  ccmpensa tors  were  not  feasible,  the  Guillemin- 
Truxal  design  procedure  was  dismissed  as  inappropriate  for  this  appli¬ 
cation. 


Entire  Eicon  structure  Assignment 

Entire  E i gens  true ture  Assignment  is  a  desl  *■-,  cue  whereby 
the  engineer  may  assign  both  the  eigenvalues  and  eigenvectors  to  a 
system.  Tiiis  technique  is  based  on  full  state  feedback  and  use  .  the 
desired  eigenvalues  and  eigenvector,  to  calculate  the  state  f.  •  s 
matrix  (Pef  10).  With  only  one  input  as  in  this  application,  t1  ■ 
selection  of  the  closed-loop  eigenvalues  completely  determines  the 
eigenvectors.  Because  of  fin's,  the  effectiveness  of  tins  technique 
was  also  questionable. 
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Presently,  the  RPV  and  the  F-15  are  different  order  systems  due  to 
tra  actuator  repp..- sen  tat  ions  of  each.  Thus,  the  use  of  this  technique 
will  he  difficult.  Tu  simplify  the  process,  it  was  decided  to  employ 
a  reduced -order  v.-rv  >  model  for  the  RPV.  The  criteria  used  for  the 
selection  of  the  r educed  order  model  v/as  time  response  to  a  pulse  input 
The  output  of  the  three  models  listed  below  was  plotted  against  the 
RPV's  second-order  model  (Fig  22)  (Eq  (6)). 
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Reduced-Order  Representation 
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Figure  22.  Reduced  Order  Servo  Appro,*, inn tion  for  the  RPV 
!odel  1  has  the  same  2  percent  setting  time  as  the  second-order  model . 


The  other  model,,  both  of  which  had  slower  responses,  were  ex ami  nod 


because  they  would  be  more  likely  to  produce  a  better  approxi  .ation  to  th 
second-order  model.  It  was  decided  tint  model  2  provided  the  most  ac¬ 
curate  approximation  to  the  second-order  model,  so  it  was  used  throughout 
the  rest  of  this  application. 

The  F-15  eigenvalues  were  assigned  as  the  closed-loop  eigenvalues 
of  the  RPV  (Ref  3).  As  expected,  the  associated  eigenvectors  did  not 
correspond  to  the  F-15  eigenvectors  (see  Table  XV).  Thus  the  Entire 
Hi  gens  true ture  Assignment  design  technique  was  also  dismissed  as  inap¬ 
propriate  for  this  application. 

Extension  cf_  Observer  Theory 

With  both  the  Gui 1 lemin-Truxal  and  the  Entire  Eigenstructure  Assign¬ 
ment  techniques  yielding  unacceptable  results,  something  mere  basic  to 
control  system  design  was  needed.  Observer  theory  gave  the  necessary 
insight.  Tin?  generalized  block  diagram  (Fig  23)  for  this  technique  is 
in  terms  of  an  actual  model,  a  desired  msuel ,  and  several  feedback 
matrices  (L,  K,  and  K  * ) .  In  an  observer,  K  and  F.1  would  be  identity 
matrices  and  the  error  feedback  would  go  to  the  desired  model  instead 
of  the  actual  model,  i.e.,  the  desired  model  would  be  made  to  track 
the  actual  model .  In  model  matching  the  goal  is  for  th -  actual  model 
to  crack  the  desired  model ;  hence,  the  error  feedback  gees  to  the  actual 
model.  K  and  K1  were  added  to  the  system  to  increase  flex  Hr  1 ity.  In 
this  application,  the  RPV  would  be  the  actual  model  and  toe  F  ■  1 3  would 
be  the  chaired  model. 

The  closed- loop  sys tom  was  expressed  in  forms  of  several  col  of 
state  variables.  It  wia  found  tint  the  best  choice  of  states  wore  the 
five  desired  model  (F-15)  state  variables  (x1)  and  the  five  error  state 
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Figure  23.  Generalized  Bloc!:  Diagram  for  Model  Matching 
variables  (c),  where 

£  =  x  -  x'  (20) 

Thus,  the  closed-loop  state  equation  is 


U  -  BLKC  |  (A  -  BLKC)  -  (A1  -  BLK’C 1 ) 


BGc  +  B1" 

+  -  u  (21) 

B1 


The  servo  representations  have  been  included  in  A  and  A1.  Also, 
for  single  input  systems,  B  and  B‘  will  always  be  zero  vectors  except 
for  the  last  element.  Hence,  K'  added  no  additional  flexibility  to  the 
system  and  was  assumed  to  be  the  identity  matrix.  For  this  application 


C  and  C'  were  also  identity  matrices.  Finally,  defining  R  -  LK  pro¬ 
vided  insight  ir.to  choosing  the  feedback  matrices.  Thus,  Eq  (21) 
became 


£ 

_ 

= 

1 

A  -  BR  !  (A  -  BR)  -  (A1  -  BL) 

I _ _ _ 

£ 

x 1 
L-J 

T 

0  1  A1 

+ 


(22) 


This  form  of  the  state  equation  yields  insight  into  the  proper 
choice  of  weighting  matrices  (R,  L,  and  Gc).  The  solution  of  Eq  (26) 
is  in  terms  of  the  eigenfunctions  (vj  e'’’t).  Thus,  R  should  be  chosen 
such  that  the  domic,  orb  terms  of  the  eigenfunctions  of  (A  -  BR)  die  out 
rapidly.  Even  though  all  eigenvalues  cf  (A  -  BR )  need  to  be  negative, 
this  does  not  indicate  tnat  all  eigenvalues  need  to  be  far  in  the  left 
hand  jjlane.  Only  the  ones  associated  with  the  dominant  eigenfunctions 
need  have  large  negative  values.  [(A  -  BR)  -  (A’  -  BL ) ]  is  a  coupling 
matrix.  Ideally,  the  best  choice  for  this  matrix  would  be  a  zero 
matrix.  Unfortunately,  L  dees  not  have  enough  variables  to  accomplish 
this.  Therefore,  1,  should  be  chosen  to  minimize  the  coupling  effect. 
Finally,  it  would  bo  desirable  if  an  input  did  not  directly  affect 
the  error.  Thus,  Gc  should  bo  chosen  to  produce  a  zero  vector  for 
(BGC  -  B'). 

The  techniques  for  choosing  R  and  Gc  are  relatively  simple. 
However,  minimization  of  the  coupling  effect  represents  a  greater 
challenge.  To  simplify  the  notation  of  this  derivation,  define 
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=  A  -  BR 


b  =  (A  -  BR)  -  (A1  -  BL) 


Thus,  Eq  (22)  becomes 


u  ,  B 

_ l _ 

| 

v  ' 

L-. 

0  |  A' 

L-  J 


;  BGC  -  B' 


B' 


(23) 


Examining  the  zero  input  response  to  Eq  (23)  produced  valuab1e  insight 
into  the  importance  of  £.  The  zero  input  response  is 


(24) 


r(t) 

10 

>it  T 

‘  l(0)“ 

= 

Z 

e  Vi  r_i 

x'(t) 

i  =  l 

x'(0) 

where 


A i  ■-  ith  eigenvalue  of 
Vi  =  ith  eigenvector  of 
r-j  =  i  th  reciprocal  basis  vector  of 


a  i 

B 

l 

- 1 

.  9  ! 

A’ 

'  a  J 

_ l_ 

(5  " 

i 

0  i 
_  1 

A' 

J 

Because  of  the  upper  block  form  of. 


:  e  i 


a  (  6 

0  !  A* 


0 


;  a- 

i  _ 


]  ,  the  system  eigenvalues  are 


simply  the  eigenvalues  of  •>  and  A'.  Also,  v,-  r^  is  a  10  x  10  matrix 
whose  upper  right  5x5  block  (defined  as  ■ )  reeds  to  be  minimized  for 

all  i. 

For  the  eigenvalues  corresponding  to  m,  the  eigenvectors  and 
reciprocal  basis  vectors  are 
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(25) 


-1 

Li 

Li  = 

wi^ 

0 

wiTDCxia  I  -  A*  3-1 

where 

vj  =  i th  eigenvector  of  a 

wj  =  i "th  reciprocal  basis  vector  of  a 


For  the  eigenvalues  corresponding  to  A1,  the  eigenvectors  and 
reciprocal  basis  vectors  are 


Y-i  = 

[■'iA'  *  _  «]  ' i 1 
'"1 

| 

Li  = 

0 

£l'  J 

wi 1 

—  — 

(26) 


v/here 

=  i^h  eigenvector  of  A' 
wi 1  =  ith  reciprocal  basis  vector  of  A' 


Now  that  the  eigenvectors  and  reciprocal  basis  vectors  have  been 
determined,  y  needs  to  be  examined.  Recall  that  L  yields  control  over 
only  the  fifth  row  of  t  .  Therefore,  the  equations  were  simplified  to 
isolate  these  terms.  For  the  eigenvalues  corresponding  to  u 


r  . 


V;  r 


ir  i 


.T 


[wiT  ;  v.'jT  ['■  iuI  -  A']'1] 


or 


y  =  :,i  wiT  Y  [Aj^1  -  A ■  ]- 1  #  (?/) 

Defining  T  ~  [  1  j  I  -  A ' ] “ ^  and  expanding  Fq  (27)  for  the  r-c  component 

of  Y 


/0 


<  rc 


k"l  '■ -ir 


i  r  wi  k 


5 

j=* 


-kj  't-jcj 


-  kIi  fir  wik 

or,  in  matrix  form 


5 

ii'l  ‘kj  i'jc 


+  f  i  rw  i  5 


“1  ^  ?jci 


Yrc  _  Lfirwi5>lc  •  •  •  firwi5>5cJ 


,:-i  'irwik  !j  =  i 


For  the  eigenvalues  corresponding  to  A1 


•51 


f‘55 


Vfti 


.T  _ 


[ A i a  1  I  ■  <<]  ^  f§i ' 


r-.i 


[0T!  Wi'Tj 


(23) 


or 


Y  =  C>.iA*  I  -  a]'1:!  ii'wi'T  (29) 

Defining  i£  =  [ > -j a '  1  -  &]_1  and  expanding  Eq  (29)  for  the  r-c  component 
Of  Y 


5  5 


Yrc  ~ 

kM  [  j  =  l 

'krjsji<]  Cik1 

wic ' 

r-. 

5  4 

k'  1  !  j  -1 

'frj;  jkJ  r i k  ' 

5 

wic '  +  k  =  ]  vr5‘ 5k 

i  k  ' w  i  c 

or,  in  matrix 

form 

1 '51 

^  1  -  r5 

>  rc  L 

ij'wic  ’  •  • 

r5  ’  i  5 ' w  \ c  1 

1  ‘  55 

5  4 

kL']  f-j^i  •  r  j  •  j  k  3  '  i  k  'wic  '  (3°) 


/I 


Aufjnontinq  Eq  (28)  with  Eq  (30) 
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re 
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sirwi5'Mc  ■  • 

•  si rwi 51 5c  ” 

or5sil'wic  ‘ 

‘  Vr5i5wic' 

*1 


4 

_k=J 

5 
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5 

'JirW'k  j  =  j  -  kj'-'jc 

4  “  ~ 

jhj  ’i'rjhjk  rH  kv/ic  ' 


51 

”55 


=  c  +  d 

where 

1  i  c  i  5 

1  s  i  c  5 

1  s  r  s  5 


(31) 


(32) 


Thus , Eq  (32)  is  a  matrix  equation  consisting  of  250  equations  in  terms 
of  5  variables. 

Minimization  of  Eq  (32)  accomplishes  the  goal  of  minimizing  the 
coupling  effect.  A  cost  function  of  one  half  the  squared  error  was 
chosen  for  this  minimization.  Thus 


J(ns)  =  \  (C  ,-_5  +  0)T  (C  ;:5  +  D) 

Setting  the  first  derivative  of  J(.:5)  equal  to  zero  and  solving  for  ,-5 
produced 

15  -  -LCrC]-1  Ct0  (33) 

Thus,  choosing  the  fifth  row  of  .■  (hence,  L)  in  accordance  with  Eq  (33) 
will  minimize  the  coupling  effect. 


For  the  model  matching  application,  22  sets  of  eigenfunctions  were 


examined  in  choosing  R  (Ref  3)  with  the  best  one  being 

R  =  [162.1  25.54  -49.12  -3.17  7.26] 

This  R  assigned  the  eigenvalues  of  a  to  :•]  =  -2,  X2  =  -4,  A3  =  -15, 

>.4  =  -17,  and  A5  =  -19  with  the  eigenfunction  corresponding  to  >.3  =  -15 
being  dominant  over  the  other  ones  by  one  order  of  magnitude.  Equation 
(33)  was  employed  resulting  in 

e5T  =  [7.89  -13.24  13.98  -4.62  -15.38] 

which  in  turn  produced 


and 


L  =  [163.41  23.33  -46.79  -3.94  2.36] 

0.99  0  0  0  0 

0  1.09  0  0  0 

K  =  0  0  1 .05  0  0 

000  0.80  0 
0  0  0  0  3.07 


In  this  application,  Gc  was  a  scalar  and  determined  to  be  3.33.  The 
above  matrices  were  substituted  into  Eg  (22)  to  produce  the  system 
state  equation  given  in  Fig  24.  Time  responses  were  produced  with 
both  the  RPV  and  the  F - 1 5  responses  on  the  same  plot  (figs  25  29) . 

These  curves  showed  reasonable  agreement,  but.  not  to  the  desired 
accuracy.  However,  with  refinement,  this  method  would  produce  “.act 
dupl i ca  t i on . 


73 


-0 . 059  0.33  -0.29  0  0  -0.046  0.25  -0.24 


VI.  Concludin']  Remark^ 

Cone  I  us  i  otis 

An  analysis  of  the  basic  RPV  provided  the  following  conclusions: 

1.  Handling  qualities  were  unsatisfied  in  the  phugeid,  dutch 
roll,  and  spiral  rrodes. 

2.  The  thrust  vectoring  unit  had  mini n:a  1  effect  cn  the  longi¬ 
tudinal  characteristics  and  r,o  effect  on  the  lateral -directional 
character!'  sties . 

3.  The  short  period  damping  ratio  and  natural  frequency,  phugoid 
damping  ratio,  and  dutch  roll  natural  frequency  varied  directly  with 
the  xcg  location.  All  other  characteristics  remained  essentially 
constant. 

The  following  conclusions  are  based  on  the  flight  control  system 
design: 

1.  Integral  cf  pitch  rate  feedback  to  the  elevator  reduced 
phugoid  oscillation. 

2.  Yaw  rate  feedback  to  the  rudder  increased  the  dutch  roll 
damping  ratio. 

3.  Yaw  rate  feedback  tn  the  ailerons  stabilized  the  spiral  mode. 

The  model  matching  application  implied  the  following  conclusions: 

1.  The  Gui 1 1 emi n-Truxal  design  technique  produced  comp  :nrators 
that  would  duplicate  only  one  state  at  a  tine.  Also,  these  compensators 
wore  impractical  to  implement. 

2.  The  Entire  Ei gens  true tjre  Assignment  technique  was  not  well 


wilted  to  si n-il e- input  systems. 


desired  duplication;  however,  i:-,L  i  i •  •  *:  4;'jri;e  of  accuracy. 

P.eco re  ’on da t ions 

Bused  on  observations  made  throughout  this  study,  the  following 
recommendations  are  posed: 

1.  Digital  Dotcom  is  a  very  generalized  algorithm  so  that  many 
different  aircraft  configurations  can  utilize  it  for  the  generation 
of  stability  derivatives.  However,  for  advanced  research  on  this  RPV 
a  more  comprehensive  algorithm  (such  as  FLEXSTA3)  might  prove  bene¬ 
ficial  for  the  model  construction. 

?..  The  cost  function  minimized  in  the  Extension  of  Observer 
Theory  was  one  half  the  squared  error.  Alternate  cost  functions 
mirj:  t  provide  better  effectiveness  in  minimizing  the  coupling  in  the 
•‘rrur  equation. 

3.  The  Extension  of  Observer  Theory  assumed  full  state  feedback 
Partial  state  feedback  is  a  more  realistic  assumption  which  could  be 
imposed  on  this  design  technique. 
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Appendix  A 
Hand  Calculations 

Several  stability  derivatives  and  all  the  control  derivatives  were 
not  computed  by  Digital  Datcom.  In  addition,  the  moments  of  inertia 
needed  to  be  generated.  Hence,  the  purpose  of  this  appendix  is  to 
derive  the  above  mentioned  terms.  All  notation  utilized  in  this 
appendix  is  defined  in  the  references  cited. 

Stability  Derivatives 

Seven  stability  derivatives  remained  undetermined:  CXu,  Cx  ,  CXa, 

& 

CXq.  CZu,  Cmu,  and  Cyr.  The  following  relationships  were  utilized  in 
their  generation  (Ref  2:19,112): 

1.  CXu  -  -2C0  -  U0  2§jl 

The  Uo  term  is  the  compressibility  effects  and  can  be  ignored. 

Thus,  for  the  TVU  attached  cases 

Cxu  =  -2Cd 
=  -2(0.05) 

=  -0.10 

and,  for  the  TVU  unattached  cases 

CXu  =  -2(0.04) 

=  -0.08 


3. 


3CQ. 

3a 


82 


Averaging  the  data  produced 


(TVU  attached) 
(TVU  unattached) 

Cx  =  0.59  -  0.37 
Aa 

=  0.22 

and,  for  the  TVU  unattached  cases 

CXa  =  0.59  -  0.35 
=  0.24 


^  =  0.37 

dOL 

f^-  =  0.35 

da 

Thus,  for  the  TVU  attached  cases 


4.  Cxq  =  0 

5.  C2„  =  -ZCl  -  U0  ^ 
iCl 

The  Uo  i—  term  also  represents  compressibility  effects  and  was  ignored. 

dU 

Thus 


cZu  =  -2CL 

=  -2  (0.59) 
=  -1.18 


6.  Cmu  =  0 

7.  Cyr  =  0 

Control  Derivatives 

Control  derivatives  are  nondimensional  derivatives  of  force  and 
moment  coefficients  with  respect  to  the  control  surfaces.  In  all,  nine 
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control  derivatives  were  needed  for  the  equations  of  motion:  Cx6e, 

CzSe’  Cn6e’  Cy&a’  Cs-<sa’  Cnsa*  Cy6r’  C*<$r’  and  Cf16r*  The  followin9 
empirical  relationships  were  used  to  compute  these  derivatives 
(Ref  15:  Chapters  10,  11,  and  12): 


Cx6e  =  0 


2.  C 


z6e  “bL<Sf  S 


=  -(3.2 0)(6.38) 
30.98 

=  -0.66 


3. 


C%e 


Ik 

c 


Hence 


■(0.66){11.72  - 

■  xcq) 

2.271 

,  %  mac 

Cm6e 

10 

-2.04 

25 

-1.94 

40 

-1.84 

4*  Cy6a=0 
5-  C£6a  =  la« ! 

=  (0. 35) (0. 38) 

=  0.13 

6.  Cn6a  =  KCL  C26a 

=  ( -0 . 20 ) (0 . 59 ) (0 . 1 3 ) 
=  -0.015 
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Hence 


7-  ("«V,Kb 


(«fi)C 


=  (5-42)(l .14) (-0.65) (1 .00) (1 .00) (6 .97) 

30.98 


=  0.90 


8.  C 


£Sr  "  Ly&r 


Zy  cos  ot  -  gv  sin  x 


„  nri  1.01  sin  0  -  (11.07  -  xcq)  sin  0 
'  u>yu  13.71 

=  0.067 

„  r  ^  ^v  cos  a  +  Zv  sin  a 

9*  Sr  "  "  Sr  - b - 

=  [(11.07  -  xcg)  cos  0  +  1.01  sin  0] 


xcq,  %  mac 

Sr 

10 

-0.42 

25 

-0.40 

40 

-0.38 

Moments  of  Inertia 

The  method  utilized  in  this  computation  consisted  of  separating 

the  RPV  into  major  components  and  calculating  their  moments  of  inertia. 

Two  transformations  and  ?  summation  determined  the  overall  moments  of 

inertia  (Ref  8:  Chapter  8).  For  accuracy  purposes,  this  derivation 

o 

was  completed  in  the  units  Ib-ir  .  .  Only  the  final  results  was 
converted  to  slug-ft^. 

Initially,  a  weight  and  balance  study  of  the  RPV’s  major  components 
had  to  be  completed,  the  results  of  which  are  in  Table  XVI.  Thus, 
the  Datcom  method  yields 
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Table  XVI.  Weight  and  Center  of  Gravity 
Locations  of  the  Major  RPV  Components 


W,  lb. 

x,  in. 

Body  under  wing 

43.65 

60.74 

Boom,  left 

8.70 

90.00 

Boom,  right 

8.70 

90.00 

Engine 

14.35 

82.00 

Fuel 

30.00 

39.00 

Horizontal  Tail 

4.00 

114.45 

Instruments 

10.00 

25.50 

Lead:  Case  1 

55.83 

18.00 

Case  2 

45.99 

18.00 

Case  3 

25.20 

18.00 

Case  4 

32.28 

18.00 

Case  5 

24.18 

18.00 

Case  6 

7.07 

18.00 

Nose 

5.00 

16.00 

TVU 

20.00 

106.30 

Vertical  Tail,  left 

5.80 

138.00 

Vertical  Tail,  right 

5.80 

138.00 

Wing 

34.00 

64.25 
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1.  Wing 


Thus 


Thus 


I  =  if  C-Ca3  +  Cb3  +  Cc2  Cb  +  Cc  Cb2  +  Cc3] 

=  [-(5.03)3  +  (29.03)3  +  (30.25)3(29.03) 

+  (30.25)(29.03)2  +  (30. 25)3] 

=  10871.69  lb-in.2 

Ift  -  Kn  t  (M.)2 
*oy  K0  1  Ww 
- 

=  0.70  10871.69  -^1^^ 

=  1492.17  lb-in.2 

T  _  Wwbw2  Ki  fCr  +  3Ct 
x0x  24  L  cr  +  Ct_ 

(34.00(164.5)2  (0.86)  ^30. 25  +  3(24) 
24  L  24  +  30.25  _ 

=  62138.60  lb-in. 2 

l0x  =  %  +  l0x 

=  1492.17  +  62138.60 

=  63630.77  lb-in. 2 

2.  Horizontal  Tail 

I  =  if  [-Ca3  +  Cb3  +  Cc2  Cb  +  Cc  Cb2  +  Cc3] 

=  [4(15.5)3] 

=  320.25  lb-in.2 


=  61.68  lb-in.2 
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C  r  +  3C  ^ 

Cr  +  Ct_ 

-  (4)(60)2(0.92)  fl 5.5  +  3(15.5)] 
24  (_  15.5  +  15.5  J 

=  1104.4  lb-in. 2 


In  = 

Ox  24 


I0Z  =  I0X  +  % 

=  1104.00  +  61.68 
=  1165.68  lb-in.2 


3.  Vertical  Tail 
_  Wy  byl2  J<5 


IOx  T8 


?  +  2  On/JLtAL  ? 

TCrv  +  Cty)^ 

=  (5.8)(26)2(1 .22)  fj  +  2(25.85)(15.6)  ~ 
18  [  (25-85  +  1 5 . 6 ) 2 

=  390.49  lb-in.2 


I  =  Tf  (-Ca3  +  Cb3  +  Cc2  Cb  +  Cc  Cb2  +  Cc3) 

=  °yP-°  [-(19.606)3  +  (25.850)3  +  (35.206)2  (25.850) 
=  2541.91  lb-in.2 


Thus 


'Oz 


*0 


(Wy  X)2 
Wv 


0.771 


2541.91 


(113.56)2 

5.8 


245.55  lb-in.2 


iQy  “  I Ox  +  lOz 


=  390.49  +  245.55 
=  636.04  lb-in.2 
4 .  Body 

Initially,  the  body  was  separated  into  four  different  components  which 
needed  to  be  consolidated.  This  was  done  with  the  following  results: 
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Weight,  lb 
xcg,  in. 
ycg,  in. 
zcg,  in. 


Cases 

Cases 

1-3 

4-6 

86.05 

66.05 

74.65 

65.06 

0 

0 

-0.79 

-2.54 

Thus,  for  cases  1-3 


Tn  _  Wf_SsJ<2 
37.68 


'3d  £B 
2£b  d 


_  (86 .05) ( 71 26 .1 1 ) (0.88) 
37.68 

=  74919.94  lb-in. 2 
2 


3(27.4)  13 5_ 
2035)  27.4 


In  =  Wf-K3 
i(Jx  4 


'  Ss 

71 


(86.05) (1 .04) 
4 


7126.1ll  2 

7TT35T. 


=  6316.29  lb-in.2 
lOz  =  lOy 

=  74919.94  lb-in.2 


For  cases  4-6 


Trt  _  (66.05) (4160.45) (0.42) 
~  37.68 
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=  17964.73  lb-in.2 

(66 .05) (0. 78) 
x  4 

=  3442.85  lb-in.2 


3  ( 1 4_.  4  7 )  .  _JT 
2(811  "  14, 


4160.45 

,rl8TT 


lOz  lOy 

-  17964.73  lb-in.2 


47 
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5.  Engine 


Assuming 


Assuming 


IOy  =  0.061  [|  Wpde2  +  Weie2  +  (Wp  -  WeHp2] 

=  0.061  [|  (14.35)  (26)2  +  (14.35)(8)2] 

=  499.82  lb-in. 2 
I0X  =  0.083  Wpde2 

=  0.083  (14.35)(26)2 
=  805.15  lb-in.2 
IOz  =  lOy 

=  499.82  lb-in.2 

6.  Fuel 

a  24  x  12  x  4  in.  tank 

lQy  =  TJ  Ux2  + 

=  (122  +  42) 

=  400  lb-in.2 

iox  =  Y2  {ly2  + 

=  y|  (242  +  122) 

=  1800  lb-in.2 

7.  Instruments 

a  3  x  4  x  6  in.  pack 

^Oy  =  yj  Ux2  +  *z2) 

I Oy  ]£  [(3)2  +  (6)2] 

=  37.50  lb-in. 2 
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r 

I0x  =  yjS-  Cfy2  +  *z2l 
=  |f  C(4)2  +  (6)2] 

=  43.33  lb-in.2 
IOZ  =  [**2  *  «y2] 

-  [(3)2  +  (4)2] 

=  20.83  lb-in.2 
8.  Lead 


Assuming  a  cubic  shape,  the  volume  of  lead  needed  and  the  length  of  a 
side  was 


Case 

W,  lb. 

Vol ,  in. ^ 

£.,  in . 

1 

55.83 

136.17 

5.14 

2 

45.99 

112.17 

4.82 

3 

25.20 

61.46 

3.95 

4 

32.28 

78.73 

4.28 

5 

24.18 

58.98 

3.89 

6 

7.07 

17.24 

2.58 

Since  a  cubic  shape  was 

assumed,  all 

the  moments  of 

inertia  will  be 

the  same. 

,0  -  ,2 

Hence 

Case 

lo,  lb-in. 2 

1 

245.83 

2 

173.08 

3 

65.53 

4 

98.55 

5 

60.98 

6 

7.84 
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Before  the  summation  of  individual  moments  or  inertia  could  be 
done,  each  moment  of  inertia  needed  to  be  transformed  into  the  RPV's 
axis  system.  The  first  transformation  was  to  a  set  of  parallel  axes 
located  at  the  RPV’s  nose.  Hence 

lx  =  £[w(y2  +  z2)  +  IOx] 

I y  =  £[W(x2  +  Z2)  +  IOy] 

Iz  =  £[W(x2  +  y2)  +  loz] 

Ixz  =  C[Wxz  f  IOxzl 

The  second  transformation  was  to  the  aircraft’s  center  of  gravity  via 

IOx (eg  =  lx  -  W  (y2  +  z2) 

!0y |Cg  =  Iy  -  W  (x2  +  z2) 

!0Z |cg  =  !z  -  W  (x2  +  y2) 

T0XZ | Cg  =  !XZ  ~  W  XZ 

Substituting  the  appropriate  values  from  Tables  XVI  and  XVII  into  the 
above  two  transformations  provided  the  RPV’s  moments  of  inertia 
utilized  in  this  study  (Table  XVIII). 
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Table  XVII.  Moments  of  Inertia  of  the 
RPV's  Major  Components 


lOx.  lb  in.2 

lOy,  lb  in.2 

l0z>  lb  in.2 

l0Xz>  lb  in.2 

Engine 

805.15 

499.82 

499.82 

0 

Fuel 

400.00 

1480.00 

1800.00 

0 

Horizontal  Tail 

1104.00 

61 .68 

1165.68 

0 

Instruments 

43.33 

37.50 

20.83 

0 

Vertical  Tail,  left 

390.49 

636.04 

245.55 

0 

Vertical  Tail,  right 

390.49 

636.04 

245.55 

O 

Wing 

62138.60 

1492.17 

63630.77 

0 

Body:  Case  1-3 

6316.29 

74919.40 

74919.40 

0 

Case  4-6 

3442.85 

17964.73 

17964.73 

0 

Lead:  Case  1 

245.83 

245.83 

245.83 

0 

Case  2 

178.08 

178.08 

178.08 

0 

Case  3 

65.53 

65.53 

65.53 

0 

Case  4 

98.55 

98.55 

98.55 

0 

Case  5 

60.98 

60.98 

60.98 

0 

Case  6 

7.84 

7.84 

7.84 

0 
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Appendix  B 

Aircraft  Equations  of  Motion 

The  aircraft  equations  of  motion  are  a  set  of  nonlinear,  coupled 
force  and  moment  equations  defined  in  inertial  space.  The  purpose 
of  this  appendix  is  to  derive  these  equations  in  their  complete  form 
and  to  simplify  them  into  a  more  useful  form. 

Definitions 

Initially,  two  fundamental  definitions  needed  to  be  made,  inertial 
velocity  and  inertial  angular  velocity.  These  two  quantities  were 
defined  in  body  fixed  coordinates  as 


V  =  U  Xb  +  V  Yb  +  W  Zb  (34) 

u>.bi  =  P  Xb  +  Q  Yb  +  R  Zb 


Applying  Newton's  second  law  in  component  form  yielded  the  three  force 
equations. 
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Moment  Equations 

The  vector  moment  equation  was  given  by  (Ref  11:286,  418). 


™  *  K  <H> 

■  K  (H)  +  ^ 


x  H 


where 


Ix  0 

-Ixz] 

P 

H  =  [Xb  Yb  Zb] 

0  Iy 

0 

1 

Q 

-IXz  0 

Iz 

R 

Expanding  Eq  (37),  substituting  into  Eq  (36),  and  simplifying 


(36) 

(37) 


EM  =  [PIX  f  (-R-PQ)Ixz  +  QR(Iz  -  Iy)]Xb 
+  [ PR ( I x  -  Iz)  +  (P2-R2)Ixz  +  Q  Iy]Yb 
+  [PQ(Iy  -  lx)  +  (QR-P)Ixz  +  R  Iz]Zb 

Recall 

IM  =  L  Xb  +  M  Yb  +  N  Zb 


Thus,  the  three  moment  equations  were 


EL  =  PIX  +  (-R  -  PQ)IXZ  +  QR( Iz  -  Iy) 
EM  =  PR(IX  -  Iz)  +  (P2  -  R2)Ixz  +Q  Iy 
zN  =  PQ( Iy  -  lx)  +  (QR  -  P)Ixz  +  R  Iz 


Assumptions 


(38) 


Equations  (35)  and  (38)  are  the  aircraft  equations  of  motion  in 


their  complete  form.  Several  assumptions  were  made  in  their  derivation 


and  are  listed  below: 


The  mass  of  the  aircraft  remains  constant. 


-  The  aircraft  is  a  rigid  body. 

-  The  earth  is  an  inertial  reference  frame. 

-  The  Xb  and  Zb  axes  form  a  plane  of  symmetry. 
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Decoupling  the  Equations  of  Motion 

Due  to  the  coupling  terms,  the  complete  form  of  the  aircraft 
equations  of  motion  was  very  difficult  to  use.  By  assuming  perturba¬ 
tions  about  straight  and  level  flight,  these  equations  decoupled  into 
a  longitudinal  and  lateral -directional  set  of  equations. 

The  longitudinal  equations  had  only  one  control  input,  a  deflection 
of  which  would  not  cause  any  P,  R,  or  V.  Hence,  Eqs  (35)  and  (38) 
simplified  to 

eFx  =  m  (U  +  QW) 

EFZ  =  m  (U  -  QU)  (39) 

EM  =  Q  Iy 

Similarly,  the  lateral -directional  equations  had  two  control 
inputs,  deflections  of  either  would  not  cause  any  Q.  Hence,  Eqs  (35) 
and  (38)  simplified  to 


EFy 

=  m  (V  +  RU  -  PW) 

EL 

=  PIx  -  RIxz 

(40) 

EN 

=  P I xz  +  Riz 

Longitudinal  Equations 

Simplification  of  Eq  (39)  was  obtained  by  assuming  stability  axes 
and  small  perturbations  about  straight  and  level  flight.  These  as¬ 
sumptions  constituted  a  linearization  about  an  arbitrary  equilibrium 
point.  Thus 

U  =  Uo  <■  u 
U  =  ii 
Q  =  q 

Q  =  q 
W  =  w 

W  =  w 
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r 


Substituting  into  Eq  (39)  and  ignoring  second  order  terms  nroduced 
the  linearized  perturbation  equations 


eFx  =  mu 

EFz  =  m(w  -q  Uo) 

EM  =  q  Iy 


(41) 


The  longitudinal  forces  and  moments  were  of  three  types:  aerodynamic, 
gravitational ,  and  thrust  forces.  Assuming  a  constant  altitude,  the 
aerodynamic  forces  consisted  of  lift  and  drag  which  were  a  function  of 
U,  a,  a,  q,  and  6e.  The  gravitational  forces  were  a  function  of  0. 
Hence,  assuming  that  the  thrust  line  was  along  Xb  and  through  the 
center  of  gravity,  the  longitudinal  forces  and  moments  could  be  ex¬ 
pressed  in  functional  form  as 


£FX  =  Fx  (U,  a,  a,  6,  q,  6e,  6j) 

IFZ  =  F2  (U,  a,  a,  0,  q,  6e)  (42) 

EM  =  M  (U,  a,  d,  q,  6e) 


Expanding  the  x-force  equation  of  Eq  (42)  in  a  first  order  Taylor  series 
expansion 


Fx  =  u  +  it*  a  +  ; 


dU 


30 


da 

3F 

0q 


3d 

♦  ,  ♦  gi  50 .  fT 


1 


(43) 


But  St  =  0  because  thrust  (i.e.  throttle  setting)  was  constant. 

d  *>y  1 

Also 


Thus 


gravity 


-mg  sin  e 


=  -mg  cos  0 
=  -mg 
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Substituting  the  above  and  Eq  (43)  into  Eq  (41) 


mu 


.  9FX  n  ,  3FX 
+  — a  +  — “■  fp 

3q  H  3<se 


(44) 


Substituting  'l  =  ^  and  the  longitudinal  stability  derivatives 
(Table  XIX)  into  Eq  (44) 


mU 


if-  '0  -  Cxu  'u  “  2%Cxd  a  -  CXa  a 


-Cxq  q  +  §  9  =  Cx5e  5e 


(45) 


2U0"*q  h T  Sq 

This  same  procedure  was  used  for  the  z-force  and  pitching  moment 
equation.  Note  the  substituteion  of  W  *  Upa  was  made  in  the  z-force 
equation.  Hence  the  remainder  of  Eq  (41)  became 


-Czu  'u  + 


mUp  _c_ 
Sq  2Uq 


CL. 


a  +  CLa  a 


_ 

2Uq  WLq  Sq 


cl 


q  8  CZ5e  6e 


(46) 


C  I 

■Cmu  'u  "  5Ti7T  Cm.  &  -  Cmn  a  +  q 


2U0 


Sqc 

Q 

2Uq  Cmq  4  =  ^m6e  5e  '^7) 

Equations  (45)  -  (47)  needed  to  be  rearranged  into  a  more  useful 
form,  i.e.,  a  state  equation.  Define 

Cl  =  #  +  onr  Cl 
Sq  2Uq 
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Table  XIX.  Definitions  of  the  Longitudinal  Stability  Derivatives 
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Solving  Eq  (46)  for  a 


•  =  ^ZU  I..  _CLg  +  1 

a  "Cl  U  ~CT  Cl 


rnllO  _c_  r 

Sq  "  2Uq  Lc> 


q  +  ^<$e  «e 


Substituting  Eq  (48)  into  Eq  (45)  and  solving  for  'u 
i  *  _  I”  Sq"  r  ,  SqT  CXa  Czul 

u  -  [  isub Cxu  W  cr^J  u 

Qy  +  ^Xa  CLal  0 

Lxa  2SU02  CT  J 


Sq_ 

mUo 


-_1 

uo 


,  Sar  CxA  mU&  cCLq  .  S^-  r  . 
L^V  "xt  _ "  ~2Uo  +  2hv0q2  CxqJ 

Sac  Cx8Cz'Se  ,  %  r  ] 

Cl  mUn  Cx5e  fie 


2mUo^  Cl  mUo 
Substituting  Eq  (48)  into  Eq  (47)  and  solving  for  q 

m 


q  = 


if-  Cmu  *  S a£.  -”#1 

ly  2UoIy 

Sqc“  r  SQ'C? 

ly  Cm-  '  2U0Iy  “Cl 


(48) 


(49) 


__2 

S'ic. 


rn»l. 


2UQly  Cl 


mUo _ c_  p. 

Tcf  '  2U0 


Sqc  C,n5e  +  Sqc 


r  2 


Cm-C  z.^e 


l  ly 


2U0Iy 


Cl 


+  Sq--  c 
2Uoiy  mq 


5e 


(50) 


Choosing 


X]  =  'u 
X2  =  a 
x3  =  0 

x4  =  q 
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Eqs  (48)  -  (50)  formed  the  desired  state  equation  as 


where 


A  =  A  x  +  B  6ecmd 


Cl 

All 

a12 

a13 

Al4 

a21 

a22 
A23 
A24 
a31 
a32 
a33 
A  34 
a41 


mUQ.  . 
Sq 


c 

2U0 


cl. 


mUQ 


Cxu 


+ 


^gc 

2mUo^ 


Cx«Czu 

Cl 


Sq  r  Sq-C*  Cx^CL„ 

mU0  Cxa  "  SDo2  Cl 


UO 


Sqr 

2niUo 

CZt| 

Cl 


Cxa  mUQ 
Cl  Sq 


2Uo 


CL  q 


0 

0 


1 

Sqc  r  Sqc2  CM^C^u 

~Iy  mu  2UQly  "  Cl 


(51) 
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a42.5£c,„  -&£  EaiEta 


[y 


2UQly  Cl 


A43  =  0 

A44  =  2U0ly 


Cm  * 

ct 

Cl 


mUo 

s:q: 


511 


mUo 


C  + 

LX6e  2mU0‘ 


r  CLq 

‘  “2Uo 
CxdCzse 

Cl 


+  C 


mq 


b21  = 


'25e 


831 

B41 


=  0 


Iy  Se  2Uq  ly 


Cl 


Lateral -Pi rectional  Equations 

As  in  the  longitudinal  equations,  a  linearization  about  an 
arbitrary  equilibrium  point  was  completed.  Thus 


U  =  Uo  +  u 

V  =  v 

V  =  y 


W  =  w 


P  =  p 
P  =  p 
R  =  r 
R  =  r 


Substituting  into  Eq  (40)  and  ignoring  second  order  terms  produced  the 
linearized  perturbation  equations 


.vFy  =  m(v  +  Uor) 

HL  *  p  Ix  -  r  IXz 
EN  =  -p  I xz  -  r  Iz 


(52) 
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The  lateral-directional  forces  and  moment s  consisted  only  of  aero¬ 
dynamic  and  gravitational  forces.  As  before,  the  aerodynamic  forces 
were  lift  and  drag;  however,  in  this  case  they  were  a  function  of  3, 
r,  p,  6a,  and  fir.  The  gravitational  forces  were  a  function  of  only  4> . 
Hence,  in  functional  form 


r.Fy  =  Fy  (3,  r,  p,  fia,  fir, 4.) 

ll  =  L  (3,  r,  p,  fia,  fir)  (53) 

IN  =  N  (b,  r,  p,  fia,  fir) 


Expanding  the  y-force  equation  of  Eq  (53)  in  a  first  order  Taylor  series 
expansion 


Recal 1 

Therefore 


Also 


r  +  p 

3p  F 


+  ify. 

3fia 


fia 


+ 

3fir 


fir  + 


3Fy 

34> 


* 


^[gravity  *  m9  sl'n  * 

if  =  ”3  cos  * 

=»  mg 


v  =•=  llo 
v  -  Uq3 

Substituting  the  above  and  Eq  (54)  into  Eq  (52) 


(54) 


mUo  ( '  +  r)  =  3  +  -  r  +  p 

•I  P  J  '  O  P 


+  fia  +  -ry^-  fir  +  mq  4 
4ca  ofir  ’ 


(55) 
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Substituting  the  lateral -directional  stability  derivatives  (Table  XX) 
into  Eq  (55) 


r  c  +  mU°  •  + 
-CY3  -  +  ^  e  + 


r.i'-'o 


S:q  ‘  2Uq  Cyr 


r  - 


mg 


2Uq  CyP  P  ~  63  +  Cy<5r 


fir 


(56) 


This  same  procedure  was  used  for  the  rolling  and  yawing  moment 
equations  with  the  following  results 


-C*S  8  '  2Uq  Clr  r  ~  Sg  r  "  2Uq  °2P  p  +  Sqb  P 


=  C£fia  5a  +  ClSr  fir 


"CnT  6  "  2U0  Cnr  r  +  5*jb  r  "  2Uq  CnP  P  *~SEfo  p 


cn5a  ^  +Cn<5r  fir 


(57) 


(58) 


As  before,  Eqs  (56)  -  (58)  needed  to  be  put  in  state  equation 
form.  Solving  Eq  (56)  for  6 


‘  *  s§  cys 


53-b_o  c  v 

2roUo2 


1 


r  +  rr^- 

uo 


+  Cv„  P  +  -  Cv_  fia  +  Cvrv.  fir 


2rnUo2  yP  P  mUo  y6a  mUo  y5r 

Solving  Eq  (57)  for  p 


(59) 


5_Pk  r 

,  +  .sgbf. 

Co  r  + 

hi.  r 

lx  -fi 

2Uo  I  x 

lx 

Sc[b2  r 

■  p  +  M 

C  o  f  A  A 

+ 

2  U  o  I  x  { 

P  P  lx 

lx 

(60) 


Substituting  Eq  (60)  into  Eq  (58)  and  solving  for  r 
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Table  XX.  Definitions  of  the  Lateral -Directional  Stability  Derivatives 


1 

C2 

'  Sqb  p  , 
L  iz  CnB 

Ixz  Sg_b  p 

TI  ‘lx‘> 

1 

[  Sqb£  c 

.  Ixz 

Sgb2 

C2 

L 2U0Iz  "r 

Iz 

2U0IX 

1 

~sak  c 

.  Ixz 

Sqb2 

C2 

2U0 I z  nP 

Iz 

2U0Ix 

1 

’Sgb  c 

+  Ixz. 

Sqb 

C2 

Iz  n6a 

Iz 

"Ix 

1 

"sgb  c 

+ 

sgb 

C2 

Iz  n>Sr 

Iz 

Ix 

■] 


6a 


(61) 


where 


C2  =  1  - 


Ixz^ 
lx  Iz 


Substituting  Eq  (61)  into  Eq  (60)  and  solving  for  p 


P  = 


Ixz' 


[  ^  TFlz  CnH  *  ^  +  C2~l¥Tz 


Sab2  Ixz 


Sqb2  Ixz  c  +  Sqb2  c  ■,  + 
_2%2  lx  Iz  LnP  2U0Ix  JP  L 


C2 


Ix 

Iz 

Ixz 

Ix 

Iz 

Ixz 

Ix 

n 

Ixz 

Ix 

17 

Cnr  + 


Ctr 


Sqb  2 

2UoIx  WJlr 
2 


1  + 


Ixz' 


C2  Ix  Iz 

_ Lxz? _ 

C2  Ix  Iz 


1  +  7T 


Ixz' 


'6r 


6r 


C2 1  x  I  z 
2 


1  +  -Liz 

i  r  rm 


C2 1  x  I  z 


6a 


6r 


(62) 
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Choosing 


X-j  -  & 


x2 


X-J  = 


x4  =  p 

Eqs  (59),  (61),  and  (62)  formed  the  desired  state  equation  as 


where 


C2  =  1 


Ixz2 

Ixlz 


A"  -5%  CH 

fl13  = 


A14  2Hjo2  C'/p 


A  21  " 


C2 


3k  c„  +  5ak  Co 

lz  +  T,  i„ 


lx 


Aoo  =  -  1  r-s.qk?.  Cn  +  Ixz  3k?_  c 
*22  C.  ,  2U0IZ  "r  Iz  2UoI x  r 


a23  =  0 


Sqb2 

2U0IZ 


+ 


Sqb2 

2u0rx 


(63) 
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A31  =  0 


A32  =  0 


A33  =  0 


A34  =  1 


a41 


_  Sgb  Ixz _  r  +  Cp 

'  C2  hdz  LH  Ix 


I  T  r  o  t 


Ixz' 


C2IxIz 


A42  =  20§l  tS 7  Cnr  +  2U^T7  Cs-r 


1  +  Ixz^ 
C2l7Iz 


A43  =  0 


A 


Sqb2 
44  "  2UqC2 


Ixz 

Ixlz 


Cnp  + 


Sgb2 

2U0Ix 


1  + 


Ixz^_ 

C2lxlz_ 


BI 1  =  SOo  Cy<5a 
BI2  =  cysr 


B21 

1 

"  C2 

'Sgb 

L  Iz 

Sa  +  X 

Sgb 

Ix 

<v 

b22 

1 

_  C2 

fsgb 
_  Iz 

Kr  * 

Sgb 

"Ix 

C?6r 

b31 

=  0 

B32  =  0 

b41  =  UI7  Cn5a  +  ^  Ci 


6  a 


1  +  _Lx_z_ 
1  C21xlz 


B42 


-  §3b  Ix z.  r  +  C 
“  C2  Ixlz  Cn6r  Ix  * 


6r 


I  + 


Xxzi_ 

C  2 1  x  I  z 
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Summary 

The  complete  form  of  the  force  and  moment  equations  can  be  found 
in  Eqs  (35)  and  (38).  The  decoupled  linearized  perturbation  equations 
can  be  found  in  Eqs  (41)  and  (52).  Lastly,  the  longitudinal  and  lateral- 
directional  state  equations  can  be  found  in  Eqs  (51)  and  (63). 
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Appendix  C 
Time  Responses 

Representative  time  responses  were  generated  for  the  basic  RPV, 
the  augmented  RPV,  and  the  nonlinear  simulation  and  are  contained  in 
this  appendix.  Each  figure  contains  a  series  of  three  time  responses, 
one  for  each  of  the  systems  mentioned  above.  This  arrangement  facili¬ 
tates  easy  comparison  between  the  different  systems.  The  longitudinal 
time  responses  are  displayed  in  Figs  30  -  35.  The  lateral  time  re¬ 
sponses  are  given  in  Figs  36  -  44.  Lastly,  the  directional  time 
responses  are  contained  in  Figs  45  -  53. 
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Figure  30.  RPV  Forward  Velocity  due  to  a  1  sec  Pulse 
Elevator  Command  of  -5  deg 
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Figure  31.  RPV  Angle  of  Attack  due  to  a  1 
Pulse  Elevator  Command  of  -5  deg 
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Figure  33.  RPV  Pitch  Rate  due  to  a  1  sec 
Pulse  Elevator  Command  of  -5  deq 
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36.  RPV  Sides!  i;-'  Angie  duo  to  a  1 
Pulse  Aileron  Cosnand  of  5  deg 


Figure  37.  RPV  Yaw  Angle  duo  to  a  1  sec 
Pulse  Aileron  Command  of  5  deg 
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Figure  38.  RPV  Yaw  Rate  due  to  a  1  sec 
Pulse  Aileron  Command  of  5  deg 
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Figure  39.  RPV  Bank  Angle  due  to  a  1  sec 
Pulse  Aileron  Command  of  5  deg 
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Figure  40.  RPV  Roll  Rate  due  to  a  1  sec 
Pulse  Aileron  Command  of  5  deg 
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Figure  42.  RPV  Aileron  Deflection  Rate  due  to  a  1  sec 
Pulse  Aileron  Command  of  5  deg 
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Figure  43.  RPV  Rudder  Deflection  Anqle  due  to  a  1  sec 
Pulse  Aileron  Command  of  5  deg 
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Figure  44.  RPV  Rudder  Deflection  Rate  due  to  a  1  sec 
Pulse  Aileron  Command  of  5  deg 


Figure  45.  RPV  Sideslip  Angle  due  to  a  1  sec 
Pulse  Rudder  Command  of  5  deg 
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Figure  46.  RPV  Yaw  Angle  due  to  a  1  sec 
Pulse  Rudder  Command  of  5  deg 
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Figure  48.  RPV  Bank  Angle  due  to  a  1  sec 
Pulse  Rudder  Command  of  5  deg 
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Figure  50.  RPV  Aileror.  Deflection  Angle  due  to  a  1 
Pulse  Rudder  Command  of  5  deg 


sec 


132 


6a  (deg/sec)  6a  (deg/sec)  6a  (deg/sec) 


fir  (deg)  fir  (deg)  fir  (dea) 


4 


Basic  RPV 


-4 


4 


Augmented  RPV 


0 


\ 

V — 


-4  4 


4 


0 


Nonlinear  Simulation 


-4  . 

0  4  8  12  16 

Time  (sec) 

Figure  52.  RPV  Rudder  Deflection  Angle  due  to  a 
Pulse  Rudder  Command  of  5  deg 
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Figure  53.  RPV  Rudder  Deflection  Rate  due  to  a  1  sec 
Pulse  Rudder  Command  of  5  deg 
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